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ABSTRACT 


Estimates  of  the  cumulative  error  associated  with  the 
third-order  Runge-Kutta  solution  of  the  two-degree  of  freedom 
equations  of  motion  are  presented.  These  estimates  constitute  a 
simple,  but  not  rigorous,  approach  to  automatic  increment  selection. 
Details  of  the  computer  program  TEA-2  which  utilize  these  equations 
are  described.  Sample  calculations  are  included.  TEA- 2  is  between 
2.  0  and  1. 7  times  faster  than  the  previous  trajectory  program  in 
addition  to  providing  greater  consistency  and  reliability. 
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INTRODUCTION 


Numerical  solutions  of  two-degree  of  freedom  trajectory 
problems  are  approximations  to  the  exact  solutions.  Intelligent 
and  efficient  employment  of  high-speed  computers  for  solving  these 
trajectory  problems  requires  that  the  errors  resulting  from  numeri¬ 
cal  methods  be  controlled.  This  report  presents  the  equations  and 
describes  the  computer  techniques  used  in  TEA-2,  a  computer  pro¬ 
gram  which  calculates  integration  step  length  automatically.  The 
integration  step  length  is  computed  such  that  the  estimated  error  in 
the  range  and/or  percentage  error  in  altitude  at  impact  compared 
to  maximum  altitude  lies  within  prescribed  limits.  This  program, 
therefore,  provides  ammunition  designers  with  the  ability  to  cal¬ 
culate  trajectories  such  that  the  cumulative  numerical  error  at  the 
impact  point  is  tolerable  and  of  pre- specified  magnitude. 

The  TEA-2  trajectory  program  requires  the  usual  input  of 
drag,  thrust,  and  mass  tables  as  well  as  initial  conditions  for  velocity, 
angle,  and  the  coordinate  position.  In  addition,  one  must  furnish 
allowable  percentage  errors  for  the  range  and/or  maximum  altitude. 
Details  concerning  input  and  output  are  included  in  Section  C  of  the 
Procedure. 

Equations  for  error  estimation  and  automatic  step  selection 
are  presented  in  Section  A.  Program  limitations  and  computer  flow 
are  described  in  Section  B.  Appendix  I  contains  a  FORTRAN  listing 
of  the  program.  Lastly,  sample  calculations  are  presented  in 
Appendix  II. 
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PROCEDURE 


A  previous  report.  Technical  Memorandum  1262,  established 
a  rigorous  mathematical  foundation  for  this  work.  Results  from 
that  report  are  utilized  in  conjunction  with  additional  approximations 
to  provide  a  set  of  equations  for  TEA-2.  Section  A  contains  this 
aspect  of  the  analysis.  The  general  features  of  TEA-2  and  detailed 
instructions  for  its  application  are  presented  in  Sections  B  tuid  C, 
respectively. 

A.  Mathematical  Basis 


Trajectory  equations  of  motion  for  a  point-mass,  two-degree 
of  freedom  system  are: 


I 


V  cos  0 

X  (t=0)  =  0 

V  sin  0 

Y  (t=0)  =  Yq 

g  sin  0  -  ^KjQd^V^/m 

V  (t=0)  = 

-g  cos  0/V 

Q  (t=0)  =  0^ 

(1) 


where  T(t),  m(t),  j’ly),  Kjj(V,  f  ),  g(y),  and  d  are  given.  They 
symbolize  thrust,  mass,  air  density,  the  drag  coefficient,  gravi¬ 
tational  acceleration,  and  particle  diameter,  respectively.  The 
quantities  V  and  0  are  the  magnitude  and  orientation  of  tbe  particle 
velocity,  respectively.  Particle  position  coordinates  are  specified 
by  X  and  y.  Figure  1  illustrates  these  quantities. 
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Figure  1 

The  Coordinate  System 


An  approximate  solution  to  the  set  of  equations  (1)  is  obtained 
by  using  a  third-order  Runge-Kutta  method.  The  approximate  solution 
is  given  by: 

where 

K,u  *  (^> 

^2U  ”  **  ^tyf  K'  ^  t 

ii  ^C,hJ 

^5U  *  ^ ^  ^*^tx  ^  ^  ^ZX i  ^  ^  ) 
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Coefficients  corresponding  to  the  Gill  modification  are  used. 
They  are: 


Cl  =  . 62653829 
C^  =  55111241 

C4  =  -.48268182 
C5  =  .  85614329 

Application  of  the  Range -Kutta  solution  method  to  equations  (1) 
introduces  truncation  error  which  is  dependent  on  the  degree  of  ap¬ 
proximation  and  rounding  error  which  follows  from  the  limitations 
of  the  digital  computer.  Errors  at  the  i-th  time  step  propagate 
forward  in  time  and  introduce  ''propagation  error.  " 

The  aforementioned  errors  were  treated  at  length  in  TM-1262. 
It  was  shown  that  a  bound  for  the  cumulative  error  can  be  obtained 
from  consideration  of  propagation,  truncation,  and  rounding  errors. 
Results  from  TM-1262  are  combined  with  assumptions  and  approxi¬ 
mations  in  this  section.  The  set  of  equations  so  derived  provide 
estimates  of  the  cumulative  error  (not  rigorous  error  bounds),  and  is 
employed  in  the  TEA- 2  program.  The  assumptions,  limitations,  and 
equations  used  in  TEA-2  follow. 

Assumptions: 


1.  The  speed  of  the  point-mass  does  not  change  by  more  than 
5%  per  time  increment. 

2.  The  orientation  of  the  velocity  vector  does  not  change  by 
more  than  2  degrees  per  time  increment. 

3.  The  difference  between  the  values  generated  by  a  single 
step  of  length  (h)  and  two  steps  of  length  ^  is  only  due  to  truncation 
and  rounding  errors. 

4.  Quantities  undergoing  arithmetic  operations  and  those  in 
storage  possess  a  relative  error  less  than  10“®. 

5.  Sine  and  cosine  routines  used  on  the  Picatinny  Arsenal 
IBM  709  have  a  relative  error  less  than  10"®. 


6 


Derivations  contained  in  TM-126Z  established  rigorous 
error  bounds.  This  report  is  not  concerned  with  rigorous  error 
bounds,  but  with  error  estimates.  Consequently,  bounding  values 
for  the  variables  used  in  TM-1262  will  be  replaced  by  local  values. 
Moreover,  estimates  for  truncation  error  based  upon  the  difference 
in  the  variables  for  a  single  step  of  length  (h)  and  two  steps  of 
length  (h/2)  will  be  employed.  Appropriately  modified  equations  • 
from  TM-1262  follow.  Estimates  of  cumulative  errors  at  the  i-th 
point  which  result  from  errors  at  the  n-th  point  are: 


=  e 


xn 


Cyl 


■yn 


where 


(c;^c^')c;  2/^ 

/  -hcKyjf)  -  c!  cj  ^ 

(./ 20J] 


max 

lATI  / 


^  K<^ 


f4 


^  ^23  ^  ^^24 


t  _ 


(5) 


(e) 


(7) 

(6) 

(9) 

Oo) 

(n) 


1 


and 


<2)  = 

^44^\/T  "  ^34^6T 

^44 

02) 

~^43  ^V7  ^  ^32  ®^er 

>  * 

'^43^VR  ^  ^93*^BR 

iC  * 

e-  sfss 
^34 

05) 

€  = 

y  ”^44^  '^‘^^34^3  \ 

04) 

= 

/  / 

05) 

A,  = 

*^33*^44  “^34^43  j  ^  ~ 

06^  J7) 

The  propagation  error  coefficients  are  given  by: 


*  cosQ  ^  ^5  *  Is//? ©I  Os) 

*^/3  *  Pc  ^  >0003 cV  09J) 

=  hZO /  0.26  -  .OOOJc  (i3>2) 

^23  =  ^-06)0.0 -.02c  ^.0OO4c^)  ^  -JO  (j9.3) 

^24  “  ^44 

=  .OJ9  -  JJOc  ^  .04 C  *  (J9.S) 

«?C^5*  Ofi.  ^.07)(j.40-.02c  ^  ,0003c^)  09^'^) 

c<^^=  J-20  4.2J  ^  .003c  0%^) 

where  ^  -  2  Kp  (20) 
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Next^one  invokes  the  assumption  that  the  difference  between 
the  values  resulting  from  a  single  step  of  length  (h)  and  two  steps  of 
length  (h/2)  is  only  due  to  truncation  and  rounding  errors.  Subtraction 
of  the  rounding  error  yields; 


A/  ■  maJ(  (h)  - y  (21- >) 

Ay  ^  max  <|yr/>;  -  0  >  .  ^2J.2) 

Ay  =  max  oy  C2I.3) 

A6  •  max  <\6(h)  - 


where  the  rotuiding  error  coefficients  may  be  approximated  by: 


=  A  9 

X  iO 

~7 

-7 

=  A9 

X  JO 

■J-  2.0  /  i>.5c]  y  jo' 

~7 

* 

4  JO 

It  follows  from  equations  (15)»  (23),  (92),  (93),  (94) 
(95)  of  FRL-TM-1262  that 

AV 

^  1 

4 

3 

A 

Ae 

^  7 

4 

% 

Thus, 

1  1 

^1  *=^55  1 

AV 

^  VAQ) 

A  of,,  yAe)4\A,\ 
7Vh^ 

(22) 


(25) 

(24) 

(25) 

(26) 
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Also, 


(27) 

C(yT 

7vh'* 

(Ztl) 

Time  increments  are  selected  such  that  the  estimated 
cumulative  errors  based  upon  local  (that  is,  n-th  step)  propagation, 
truncation,  and  rounding  errors  satisfy  specified  allowable  error 
criteria  for  the  x  and/or  y  variables  at  the  impact  point.  That  is. 


l^(h)  =  e„/(h)  -  ij,  -  O  (28) 

and/or  -  iy  =  0  (29) 

where  the  subscript  f  denotes  the  impact  point  and  the  quantities  Cy 
and  are  the  allowable  errors  in  the  x  and  y  variables  at  impact, 
respectively.  Equations  (7)  and  (8)  require  the  number  of  steps 
between  the  n-th  and  f-th  points.  Moreover,  estimates  of  the  range  and/or 
maximum  altitude  will  be  required  since  the  allowable  error  is  Specified 
as  a  certain  percentage  of  these  quantities. 


Equations  (E8)  and/or  (29)  are  solved  by  the  method  of  double 
false  position  with  a  convergence  criterion  of  h,  -  ^  ,  q£ 

» 

where  and  ^ ,  i  The  starting  value  of  the 

iteration  followed  by  3.S/io~  h2  if  or  by 

,./hc  -  if  >  O  ,  The  maximum  step  size  is  .035  V/g  which 

is  consistent  with  assumption  2.  Time  increments  less  than  .000035V/g 
are  not  permitted. 


The  following  formulas  are  employed  to  obtain  the  required 
estimates: 


where 


(f-r,)  =  M/h 

jir,  * 


C30) 

C3i) 


or 


At  =  -  in  such  hhof  y(t*) 
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where  is  the  starting  point  for  the  terminal  ballistic  phase  of  the 
trajectory.  The  time  t^^^  will  be  specified  in  the  thrust  table. 


The  terminal  ballistic  phase  (Phase  IV)  impact  range  and 
maximum  altitude  are  estimated  from  a  linear  drag  model.  They  are 
given  by: 


AUj 

A(Ui) 


L3Z) 


i  =  n 


where 


)/rnax  “ 

iup.  ^>1 

II 

/  AOi 

Sin  Qfj 

Je-fin€S 

Mu)  = 

AUi  ' 

{('On*. 

)  rci„  ^ . 

L32.1) 

iS3) 

IS3.I) 

C33.2) 


''  /*  // 
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y.  M.  isn  Qn 
9 


C33.3) 


u  ^  Uf., 


(.33.4) 


The  following  approximation  is  made  for  the  model 
trajectory: 


where 


^  =  [.  00234  Ib/ft^ 

]  X  10-™ 

m=0 

0  ft  ^ 

y 

58,  000  ft 

m=l 

58,  000  ft  < 

y 

100,  000  ft 

m=2 

100,000  ft  ^ 

y 

160,  000  ft 

m=3 

160,000  ft^ 

y 

230,  000  ft 

m=4 

230,  000  ft  « 

y 

300,  000  ft 

m=5 

300,  000  ft  < 

y 

400,  000  ft 

m=6 

400,  000  ft  ^ 

y 

600,  000  ft 

m=7 

600,000  ft  ^ 

y 

Estimates  of  the  dependent  variables  for  an  intermediate 
ballistic  phase  (Phase  II)  from  tj  to  tg  follow: 

=  Vj  -  ^  (.34^ 


i7  =  0.5  (Vj  LB4.I) 

6),  =  O.SCO^^Qj)  iSSj) 

*  O.Si/,cosG,(i^’-ij)  L3e) 

=  yj  i-  0.8y,  smQfd^- ij)  (^7) 

The  velocity  change  for  a  thrust  phase  from  t}^  to  tj^  has  been 
calculated  from  linear  drag  approximation.  Estimates  for  thrust  phases 
are: 

<?C  ®  ^stnQ^  (38) 

/S  =  pf<ocl^yh 


=  r(u)  -  m(iK) 

’  r  /o(* 

y- 

P  p’-th 

. 

/  m(iK)  j  -  rdt,) 

/S-m  ^  p  . 


mdx) 
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where 


/a' 


T  -  Th  * 

ro  ^ 

0.S  ()i  /  l//f) 

in) 

-  $C059^(iK-ih)  / 

di) 

0. 5  ) 

(.43.1) 

/  0.6  )/2  cos 

(.44) 

14S) 

B.  Computer  Technique 


A  number  of  refinements  have  been  incorporated  in  TEA-2  in 
order  to  generate  reasonable  error  estimates.  These  refinements 
include: 


1.  Conservative  range  estimates. 

2.  Insuring  that  the  truncation  error  coefficients  do  not 
decrease  by  a  factor  of  more  than  .  8  from  the  last  time  they  were 
computed  unless  a  phase  change  occurs. 

3.  Averaging  the  estimated  error  at  the  (n+N)-th  step  such 
that  it  is  one-half  the  sum  of  estimates  based  upon  local  values  at  the 
n-th  and  (n+N)-th  steps  provided  that  N  steps  have  been  taken  since  the 
last  step  selection  and  that  a  phase  change  does  not  occur  at  n+N. 

4.  The  acceptable  percentage  error  is  doubled  if  the 
criteria  cannot  be  satisfied. 


The  TEA-2  trajectory  program  may  be  utilized  to  compute 
two-dimensional,  point  mass  trajectories  for  one  and  two  stage 
rockets  and/or  ballistic  projectiles  and  consists  of  the  following 
phases: 


Phase  I 
Phase  II 
Phase  III 
Phase  IV 


Acceleration  of  Booster  and  Main  Stage 
Coasting  of  Main  Stage 
Acceleration  of  Main  Stage 
Free-flight  of  Main  Stage 


Any  of  the  above  phases  may  be  excluded  in  the  computations; 
for  ballistic  trajectory  computations,  only  Phase  IV  is  required. 


The  following  are  program  conditions: 


1. 

2. 

3. 

4. 

5. 

6. 

7. 

8. 


9. 


Control  systems  are  not  considered 
A  flat  earth  model  is  assumed 
Winds  are  not  considered 
ARDC  Standard  Atmosphere,  1959 
Constant  thrust  or  thrust  variable  with  time 
Thrust  modification  for  atmospheric  pressure  changes 
Constant  linear  weight  change  or  weight  variable  with 
time  for  rocket  boosted  phases 
kjQ  or  Cju  drag  coefficients  may  be  used 

Form  factor  to  alter  drag  coefficient 


Thrust  values  for  the  acceleration  phases  (I  &  III)  are  obtained 
in  the  computations  by  linear  interpolation  of  a  thrust  (lb)  versus 
time  (seconds)  table  when  table  values  are  presented;  however,  a 
constant  thrust  force  may  be  presented  instead  of  tables. 


A  constant  weight  change  (burning  rate,  lb/ sec)  or  a  table  of 
weight  (lb)  versus  time  (seconds)  may  be  introduced  into  the  computa¬ 
tions  for  rocket  boosted  missiles.  Using  the  table,  weight  is  obtained 
in  the  program  by  linear  interpolation.  When  tables  are  not  presented, 
the  weight  is  decreased  during  thrust  phases  by  a  factor  of  integration 
increment  (seconds)  multiplied  by  the  burning  rate  for  each  step  of 
integration. 
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Tables  of  drag  coefficients  (Cj)  or  kj^)  versus  mach  numbers 
may  be  presented  for  each  acceleration  phase  to  be  computed  and  for 
the  free-flight  phase.  If  only  Phase  IV  (i.  e.  ,  ballistic  trajectory) 
is  to  be  computed,  then  only  the  drag  coefficient  table  for  free-flight 
is  required. 

A  constant  drag  coefficient  may  be  presented  for  the  entire 
trajectory.  (See  description  for  Card  #1.  ) 

A  generalized  and  a  detailed  flow  chart  are  presented  in  order 
to  depict  the  program  logic.  The  corresponding  FORTRAN  listing 
can  be  found  in  Appendix  I.  Reference  should  be  made  to  the  flow 
charts  prior  to  any  modification  due  to  the  complexity  of  the  branching 
operations  in  the  body  of  the  program.  This  program  requires  7024 
words  of  core  storage. 
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C.  Computer  Usage 


The  procedures  for  preparing  input  data  to  trajectory  program 
TEA-2  are  discussed  in  this  section,  and  examples  of  typical  cases 
are  presented  and  described. 

Input  data  is  punched  on  cards  which  are  loaded  with  the  program 
deck  onto  tape.  Data  sheets  which  may  be  used  for  presenting  required 
data  for  key-punching  are  illustrated  in  the  descriptions  of  sample 
trajectory  cases  to  be  computed.  The  computer  program  reads  input 
data  from  Tape  A2  which  is  the  standard  for  IBM  FORTRAN  Monitor 
Systems. 

Descriptions  of  input  requirements  will  consist  of  two  parts: 

(1)  Rocket  trajectories  and,  (2)  Ballistic  trajectories. 
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I .  Rocket  Trajectory  Input 


Card  #1  contains  nine  one-digit  and  three  two-digit  control 
parameters  from  columns  1  through  15.  Columns  16  to  20  are  for 
a  decimal  input  number.  These  parameters  are  to  be  presented 
across  card  #1  as  follows: 

Card  #1 _ NEXT(l),  NEXT(2),  NEXT(3),  NEXT(4),  KDCON, 

NTHRST,  KD,  M,  NXY,  JN,  IMJR,  NS,  GAMMA 
Format  (911,  312,  F5.  0) 

where: 

NEXT(l)  ....  equal  to  0  for  not  computing  Phase  I 
equal  to  1  for  computing  Phase  I 

NEXT(2) ....  equal  to  0  for  not  computing  Phase  II 
equal  to  1  for  computing  Phase  II 

NEXT (3)  . . ..  equal  to  0  for  not  computing  Phase  III 
equal  to  1  for  computing  Phase  III 

NEXT (4). . . .  equal  to  0  for  not  computing  Phase  IV 
equal  to  1  for  computing  Phase  IV 

KDCON  ....  equal  to  1  for  variable  drag  coefficients 
(use  tables) 

equal  to  2  for  constant  drag  coefficient 

NTHRST. . . .  equal  to  1  for  using  thrust  tables 

equal  to  2  for  not  using  thrust  tables 


KD . equal  to  1  for  using  kj^  coefficients 

equal  to  2  for  using  Cjj  coefficients 

M . equal  to  1  for  using  weight  table 


equal  to  2  for  not  using  weight  table 
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NXY 


equal  to  1  for  step  selection  based  on  range 
error  criterion 

equal  to  2  for  step  selection  based  on 
maximum  altitude 

equal  to  3  for  step  selection  based  on  total 
range  and  maximum  altitude 

JN .  number  of  trajectories  to  be  computed  for  the 

same  conditions  of  error  limit,  thrust  and 
drag  coefficients  (a  two  digit  integer) 


NJR  ......  trajectory  identification  number  for  first 

trajectory  to  be  computed  (  a  two  digit  int.eger) 

NS .  number  of  integration  steps  to  be  computed 


for  error  estimations  (a  two  digit  integer;  if 
left  blank  on  card,  set  equal  to  ten  in  program) 

GAMMA  ...  percent  error  limit  for  trajectory  computations 
(if  left  blank  on  input  card,  GAMMA  is  set  equal 
to  5.  0%  in  program) 

Cards  #2  to  #7  contain  the  drag  coefficient  table  for  Phase  I.  Cards 
#2  to  #4  consist  of  the  table  mach  numbers,  and  cards  #5  to  #7  contain 
the  corresponding  drag  coefficients. 

Cards  #8  to  #13  contain  drag  coefficient  tables  for  free-flight  phases 
(Phases  II  and  IV).  Cards  #8  to  #i0  are  for  the  table  mach  numbers,  and 
cards  #11  to  #13  are  for  the  corresponding  drag  coefficients. 

Cards  #14  to  #19  contain  the  drag  coefficient  for  Phase  III.  Cards 
#14  to  #16  are  for  the  table  mach  numbers,  and  cards  #17  to  #19  are  for 
the  corresponding  drag  coefficients. 

Either  kj^  or  Cj^  coefficients  may  be  used  in  the  tables;  however, 
values  that  are  to  be  used  must  be  specified  by  the  control  parameter 
KD  in  card  #1.  All  of  the  coefficient  tables  may  be  eliminated  as  input 
when  in  card  #1  KDCON  =  2, 


Cards  #20  to  #25  contain  the  thrust  table  for  Phase  I,  This 
table  is  not  required  when  a  constant  thrust  force  is  to  be  used;  and 
for  this  case,  the  control  parameter  NTHRST  in  card  #1  must  equal 
2,  and  a  thrust  value  must  be  presented  in  card  #46,  Cards  #20  to 
#22  are  for  the  table  time  values,  and  cards  #23  to  #25  are  for  the 
corresponding  table  thrust  values. 

Cards  #26  to  #31  contain  the  weight  table  for  Phase  I  and 
must  be  presented  only  if  in  card  #1  M=l.  Cards  #26  to  #28  are  for 
table  time  values,  and  cards  #29  to  #31  are  the  corresponding  table 
values  of  weight. 

Cards  #32  to  #37  contain  the  thrust  table  for  Phase  III.  These 
cards  are  not  required  when  a  constant  thrust  value  is  to  be  used;  and 
for  this  case,  NTHRST  must  equal  2  in  card  #1  and  a  constant  thrust 
value  must  be  presented  in  card  #46.  Cards  #32  to  #34  are  for  table 
time  values,  and  cards  #35  to  #37  are  for  the  corresponding  table 
thrust  values. 

Cards  #38  to  #43  contain  the  weight  table  for  Phase  III  and  must 
be  presented  when  constant  burning  rates  are  not  used  (in  card  #1  M=l). 
Cards  #38  to  #40  are  for  table  time  values,  and  cards  #41  to  #43  are  the 
corresponding  table  values  of  weight. 

Data  field  widths  for  all  input  valves  are  nine  columns,  and 
eight  values  may  be  introduced  on  one  card.  Tables  are  required  only 
for  phases  that  are  to  be  computed.  Each  table  that  is  required  by  the 
program  will  consist  of  six  input  cards;  therefore,  if  data  is  not  available 
for  completing  the  maximum  table  requirements,  blank  cards  must  be 
used  to  supply  the  proper  number  of  cards. 

The  last  three  cards  required  for  input  (cards  #44,  45,  and  46) 
will  be  listed  along  with  input  symbols,  and  these  cards  are  required 
for  each  rocket  trajectory  to  be  computed.  Using  the  same  error 
criteria,  thrust,  weight,  and  drag  coefficient  tables,  "JN"  trajectories 
may  be  computed  by  varying  initial  conditions,  configuration,  etc. ,  in 
cards  #44,  45  and  46. 
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Card  #44 


THDEGO,  VO,  WGT,  DIA,  XO,  YO,  TO,  YF 


where: 

THDEGO  . . .  initial  quadrant  elevation,  degrees 


VO . initial  velocity,  feet  per  second 

WGT .  total  weight,  pounds 

DIA .  missile  diameter,  inches 

XO . initial  range,  feet 

YO . initial  altitude,  feet 

TO . time  at  start  of  trajectory,  seconds 

YF . terminal  altitude,  feet 

Card  #45 . THRST,  DOTM,  AREA,  BWGT,  THR2,  DOTM2, 


AREA2,  PRESO 

WHERE: 


THRST . booster  thrust,  pounds 

DOTM . booster  burning  rate,  lb  per  second 

AREA .  nozzle  throat  area  of  booster,  sq.  in. 

BWGT . booster  weight  (w.  o.  propellant),  pounds 

THR2 . main  stage  thrust,  pounds 


DOTM2  ....  main  stage  burning  rate,  lb  per  second 
AREA2  ....  nozzle  throat  area  of  main  stage,  sq.  in. 
PRESO  ....  static  test  atmospheric  pressure 

AREA,  AREA2  and  PRESO  are  to  be  presented  for  input  only 
when  the  thrust  is  to  be  modified  in  the  computations  for  atmospheric 
pressure  changes.  For  this  correction  not  to  be  performed,  AREA 
and  AREA2  must  be  made  equal  to  zero  or  left  blank  on  the  input 
card. 
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Card  #46 


TIM(l),  TIM(2),  TIM(3),  DWRT(l),  DWRT(2), 
DWRT(4),  CKD,  FACT 


where: 


TIM(l) .  end  of  Phase  I,  seconds 

TIM(2)  . .  end  of  Phase  II,  seconds 

TIM(3) .  end  of  Phase  III,  seconds 

DWRT(l) .  output  spacing  for  Phases  I  and  III,  seconds 

DWRT(2) .  output  spacing  for  Phase  II,  seconds 

DWRT(4) .  output  spacing  for  Phase  IV,  seconds 

CKD .  constant  drag  coefficient 

FACT . form  factor  (if  left  blank  on  input  card, 

form  factor  is  set  equal  to  one) 


ALL  DATA  PRESENTED  ON  CARDS  #2  TO  #46  MUST  HAVE 
DECIMAL  POINTS  AND  HAVE  FIELD  WIDTHS  OF  NINE  COLUMNS. 


Rocket  Trajectory  Sample  Input  Case 

The  rocket  sample  case  is  for  a  one  stage  rocket  which  would 
require  computing  Phases  I  and  IV  of  the  program. 

The  following  conditions  represent  the  trajectory  input  requirements: 

1.  Cj)  drag  coefficient  tables 

2.  constant  burning  rate  (5  Ibs/sec) 

3.  constant  thrust  (100  lbs) 

4.  step  selection  based  on  range  error  criteria 

5.  10  integration  steps  per  estimate 

6.  quadrant  elevatioii  =  45.  0° 

7.  initial  velocity  =  1000.  0  fps 

8.  weight  =  100.  0  lbs 

9.  diameter  =  11.0  inches 

10.  burning  time  =  4.0  seconds 

11.  output  spacing  every  At 

12.  form  factor  =  1.0 

13.  5.  0%  error  limit 


For  the  above  conditions,  the  following  data  must  be  presented 
in  card  #1 : 


NEXT(l) 

= 

1,  compute  Phase  I 

NEXT(2) 

= 

0,  do  not  compute  Phase  II 

NEXT(3) 

= 

0,  do  not  compute  Phase  III 

NEXT(4) 

= 

1 ,  compute  Phase  IV 

KDCON 

= 

1,  variable  drag  coefficients  presented 

NTHRST 

= 

2,  constant  thrust 

KD 

= 

2,  using  Cj)  drag  coefficients 

M 

= 

2,  not  presenting  weight  tables 

(constant  burning  rate) 

NXY 

= 

1,  step  selection  based  on  range  error 

criteria 

JN 

= 

01,  compute  one  trajectory 

NJR 

= 

01,  trajectory  numbered  1 

NS  (blank) 

,  ten  integration  steps  per  estimate 

GAMMA  (blank) 

,  5.  0%  acceptable  error 

An  illustration  of  the  presentation  of  data  on  card  #1  may  be  seen 
on  page  33. 

The  drag  coefficient  tables  for  Phase  I  and  Phase  IV  as  presented 
for  key-punching  are  shown  on  the  table  data  sheet  on  page  32,  These 
tables  are  punched  onto  cards  #2  to  #13,  and  the  punched  cards  are 
illustrated  on  page  33. 

Cards  #14  to  #19  are  not  required  for  input  as  Phase  III  is  not 
to  be  computed. 

Cards  #20  to  #43  are  not  required  for  input  as  thrust  and  weight 
tables  are  not  to  be  presented. 

The  following  data  must  be  punched  on  card  #44: 

THDEGO  (45.  0),  quadrant  elevation,  degrees 
VO  (1000.  0),  initial  velocity,  fps 
WGT  (100.  0),  missile  weight,  lbs 
DIA(ll.O),  missile  diameter,  inches 

Values  for  XO,  YO,  TO,  and  YF  do  not  have  to  be  punched  on 
the  card  for  this  case  as  they  equal  zero. 
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The  following  data  must  be  punched  on  card  #45: 

THRST  (100.0),  constant  thrust  force,  lbs 

DOTM  (5.0),  burning  rate,  Ibs/sec 

Values  for  AREA,  BWGT,  THR2,  DOTM2,  AREA2,  and  PRESO 
do  not  have  to  be  punched  on  the  card  as  this  data  is  not  required  for 
the  given  rocket. 

In  card  #46,  only  the  following  need  be  punched: 

TIM(l)  (4.0),  end  of  Phase  I,  seconds 

TIM(2)  and  TIM(3)  are  not  required  as  Phases  II  and  III  are  not 
being  computed. 

DWRT(l),  DWRT(2),  and  DWRT(4)  are  left  blank  as  output  is  to 
be  printed  for  every  integration  increment,  At. 

CKD  is  left  blank  as  table  values  of  drag  coefficients  are  presented. 

FACT  is  left  blank  as  a  form  factor  equal  to  1.  0  is  required; 
however,  a  value  of  1.0  may  be  punched  for  FACT. 

Data  presented  for  key-punching  onto  cards  #44,  #45  and  #46 
may  be  seen  on  the  input  sheet  on  page  3  1.  An  illustration  of  this  data 
on  cards  may  be  seen  on  page  33. 
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INPUT  FOR  TRAJECTORY  PROGRAM  TEA-2 


For  Ballistic  and  Rocket  Trajecrory 


QUADRANT  ELEVATION,  degrees  "^-5. 

INITIAL  VELOCITY,  ft/ sec  /OOC» 

TOTAL  WEIGHT,  lbs  /GG. 

DIAMETER,  inches  JAQ 

INITIAL  RANGE,  ft  _ 

INITIAL  ALTITUDE,  ft  _ 

INITIAL  TIME,  sec 

TERMINAL  ALTITUDE,  ft  . .. 


For  Rocket  Trajectory 

BOOSTER  THRUST,  lbs  /OO. 

BOOSTER  BURNING  RATE,  lb/ sec  S.O 

BOOSTER  NOZZLE  AREA,  sq.  inches  _ 

BOOSTER  WEIGHT  (EMPTY),  lbs  _ 

MAIN  STAGE  THRUST,  lbs  _ 

MAIN  STAGE  BURNING  RATE,  lbs/ sec  _ 

MAIN  STAGE  NOZZLE  AREA,  sq.  inches  _ 

STATIC  TEST  ATMOS.  PRESSURE,  lb  sec^ft"^  . 


For  Rocket  Trajectory 

TIME  AT  END  OF  PHASE  I,  sec 

TIME  AT  END  OF  PHASE  II,  sec  _ 

TIME  AT  END  OF  PHASE  III,  sec  _ _ 

OUTPUT  SPACING  FOR  PHASES  I  &  UI,  sec  _ 

OUTPUT  SPACING  FOR  PHASE  II,  sec  _ 

OUTPUT  SPACING  FOR  PHASE  IV,  sec  _ 

CONSTANT  DRAG  COEFFICIENT  _ 

FORM  FACTOR  _ 


For  Ballistic  Trajectory 

OUTPUT  SPACING,  sec 
FORM  FACTOR 

CONSTANT  DRAG  COEFFICIENT 
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Name 


Table  Data  Sheet  for  TEA- 2  Trajectory  Program 

Date 


2. 


Ballistic  Trajectory  Input 


Card  #1  contains  the  same  values  as  in  card  #1  described  for 
rocket  trajectories  on  page  24. 

Card  #2  to  #7  contain  the  drag  coefficient  table.  Card  #2  to  #4 
are  for  the  table  mach  numbers,  and  cards  #5  to  #7  are  for  the  cor¬ 
responding  drag  coefficients.  This  table  is  not  required  when 
KDCON  =  2  in  card  #1.  For  this  case,  the  constant  drag  coefficient 
must  be  presented  in  card  #9. 

Either  kjj  or  C^  drag  coefficients  may  be  used,  and  values 
employed  must  be  specified  by  the  control  parameter  KD  in  card  #1. 

Card  #8 .  THDEGO,  VO,  WGT,  DIA,  XO,  YO,  TO,  YF 

where: 

THDEGO .  initial  quadrant  elevation,  degrees 


VO .  initial  velocity,  feet  per  second 

WGT . total  weight,  pounds 

DIA . missile  diameter,  inches 

XO . initial  range,  feet 

YO . initial  altitude,  feet 

TO . time  at  start  of  trajectory,  seconds 

YF . terminal  altitude,  feet 

Card  #9 . DWRT(4),  FACT,  CKD 

where: 

DWRT(4) .  output  spacing,  seconds 

FACT . form  factor 

CKD .  constant  drag  coefficient 


ALL  DATA  PRESENTED  ON  CARDS  #2  TO  #9  MUST  HAVE 
DECIMAL  POINTS  AND  HAVE  FIELD  WIDTHS  OF  NINE  COLUMNS. 
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Ballistic  Trajectory  Sample  Input  Case 


The  following  conditions  are  for  the  sample  ballistic  trajectory: 

1.  Cq  drag  coefficient  table 

2.  initial  velocity  =  1000.  fps 

3.  quadrant  elevation  =  45.0° 

4.  weight  =  100.0  lbs 

5.  diameter  =  11.0  inches 

6.  output  spacing  every  ^  t 

7.  form  factor  =  1.0 

8.  ten  integration  steps  per  estimate 

9.  0.  5%  error  limit 

For  the  above  conditions,  the  following  data  must  be  presented 
in  card  #1: 


NEXT(l) 

=  0. 

NEXT(2) 

=  0, 

NEXT(3) 

=  0, 

NEXT(4) 

=  1. 

KDCON 

=  1. 

NTHRST 

=  2, 

KD 

=  2, 

M 

=  2. 

NXY 

=  2, 

JN 

=  01 

NJR 

=  01 

NS 

=  10 

GAMMA 

=  0. 

do  not  compute  Phase  I 
do  not  compute  Phase  II 
do  not  compute  Phase  III 
compute  Phase  IV 

variable  drag  coefficients  presented 
no  thrust 

using  Cj)  drag  coefficients 

not  presenting  weight  tables 

step  selection  based  on  altitude  error 

criteria 

compute  one  trajectory 


percent  acceptable  error 


An  illustration  of  the  presentation  of  data  on  card  #1  may  be 
seen  on  page  37  . 


The  drag  coefficient  table  for  the  ballistic  trajectory  as  presented 
for  key-punching  is  shown  on  the  table  data  sheet  on  page  38  .  This 
table  is  punched  onto  cards  #2  to  #7,  and  the  punched  cards  are  illustrated 
on  page  39  • 
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The  following  data  must  be  punched  on  card  #8; 

THDEGO  {45.  0),  quadrant  elevation,  degrees 
VO  (1000.0),  initial  velocity,  fps 
WGT  (100.0),  missile  weight,  lbs 
DIA  (11.0),  missile  diameter,  inches 

Values  of  XO,  YO,  TO,  and  YF  do  not  have  to  be  punched  on 
the  card  for  this  case  as  they  equal  zero. 

Card  #9  may  be  left  blank.  No  value  is  required  for  CKD  as  a 
table  of  drag  coefficients  is  to  be  presented.  No  value  is  required 
for  DWRT(4)  as  output  is  required  every  integration  step,  Ut.  No 
value  is  required  for  FACT  as  the  form  factor  must  equal  one. 

Data  presented  for  key-punching  onto  cards  #8  and  #9  may  be 
seen  on  the  input  sheet  on  page  37  .  An  illustration  of  this  data  on 
cards  may  be  seen  on  page  38  . 
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INPUT  FOR  TRAJECTORY  PROGRAM  TEA-2 


For  Ballistic  and  Rocket  Trajectory 


QUADRANT  ELEVATION,  degrees 
INITIAL  VELOCITY,  ft/ sec 
TOTAL  WEIGHT,  lbs 
DIAMETER,  inches 
INITIAL  RANGE,  ft 
INITIAL  ALTITUDE,  ft 
INITIAL  TIME,  sec 
TERMINAL  ALTITUDE,  ft 


For  Rocket  Trajectory 


BOOSTER  THRUST,  lbs 
BOOSTER  BURNING  RATE,  lb/ sec 
BOOSTER  NOZZLE  AREA,  sq.  inches 
BOOSTER  WEIGHT  (EMPTY),  lbs 
MAIN  STAGE  THRUST,  lbs 
MAIN  STAGE  BURNING  RATE,  lbs/ sec 
MAIN  STAGE  NOZZLE  AREA,  sq.  inches 
STATIC  TEST  ATMOS.  PRESSURE,  lb  sec^ff'* 

For  Rocket  Trajectory 

TIME  AT  END  OF  PHASE  I,  sec 
TIME  AT  END  OF  PHASE  II,  sec 
TIME  AT  END  OF  PHASE  III,  sec 
OUTPUT  SPACING  FOR  PHASES  I  &  III,  sec 
OUTPUT  SPACING  FOR  PHASE  U,  sec 
OUTPUT  SPACING  FOR  PHASE  IV,  sec 
CONSTANT  DRAG  COEFFICIENT 
FORM  FACTOR 

For  Ballistic  Trajectory 


-^c>^ 

_ 


OUTPUT  SPACING,  sec 
FORM  FACTOR 

CONSTANT  DRAG  COEFFICIENT 
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Name 


Table  Data  Sheet  for  TEA-^  Trajectory  Program 
_______  Date 


BALLISTIC  TRAJECTORY  INPUT 


RESULTS  AND  DISCUSSION 


The  TEA-2  digital  computer  program  for  the  IBM  709  has 
been  tested  on  a  substantial  number  of  trajectories.  This  program 
has  always  selected  reasonable  time  steps.  Time  steps  for  TEA-E 
are  dependent  upon  the  accuracy  specified  by  the  user  as  well  as 
the  problem  characteristics.  That  is,  problem  initial  conditions, 
body  drag  characteristics,  and  vehicle  thrust  characteristics  are 
important  parameters  and  exert  a  significant  influence  on  the  time 
increments  used  by  TEA-2. 

Analysis  of  output  for  typical  ballistic  flights  has  shown  that 
the  estimated  range  error  is  somewhat  more  than  10  times  too  con¬ 
servative.  That  is,  an  estimated  range  error  of  1%  by  TEA-2  would 
correspond  to  an  actual  range  error  less  than  0.  1%. 

Four  typical  trajectories  were  computed  in  order  to  establish 
timing  estimates.  The  program  limitation  to  a  maximum  angular 
change  per  time  step  of  two  degrees  places  an  upper  limit  on  compu- 
tational  speed.  This  limitation  explains  the  small  difference  in 
computation  time  between  the  5%  and  the  0.  5%  ballistic  trajectories 
presented  in  Table  1. 

The  first  ballistic  trajectory  corresponds  to  the  sample  input, 
and  sample  output  presented  in  Section  C  and  Appendix  II,  re¬ 
spectively.  A  similar  trajectory  and  the  time  steps  required  for  less 
than  0.  5%  range  error  are  presented  in  Figure  2.  Actual  range  error 
is  probably  much  less  than  0.  5%  in  this  case.  Figure  2  illustrates 
that  maximum  time  steps  were  taken  from  10.  483  seconds  to  33.  837 
seconds. 

Figure  3  illustrates  a  typical  rocket  trajectory  which  is 
similar  to  the  ballistic  case.  Initial  step  size  is  much  smaller  in 
this  case.  The  average  step  size  is  25%  less  than  for  the  ballistic 
case. 
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Table  1  lists  two  typical  ballistic  trajectories  and  two  typical 
rocket  trajectories.  Both  5%  and  0.  5%  range  accuracies  were  re¬ 
quested  for  these  trajectorie  s.  These  problems  were  run  both  on 
TEA-2  and  on  the  "standard"  two-degree  of  freedom  program  which 
utilizes  a  fixed  time  increment  which  is  an  input  quantity.  This 
input  time  increment  is  based  on  the  user's  experience.  The  values 
reported  in  the  table  are  typical  of  those  presently  used.  Under  the 
assumption  that  the  four  cases  presented  with  flight  history  con¬ 
stitute  the  "average"  problems,  the  use  of  TEA-2  constitutes  a  time 
savings  of  2.  1  and  1. 7  for  the  5%  and  0.  5%  accuracy  specifications, 
respectively.  This  will  result  in  a  cost  reduction  of  35,  900  dollars 
per  year  or  28,  500  dollars  per  year  for  the  respective  required 
accuracies.  Cost  savings  are  based  on  current  usage  of  19  hours 
per  month  at  a  cost  of  300  dollars  per  hour.  The  5%  requested  ac¬ 
curacy  should  be  the  most  widely  used  case. 

Computation  time  is  approximately  2,  5  to  3.  5  times  faster 
than  the  "standard"  program  for  the  usual  ballistic  and  rocket 
trajectories.  However,  TEA-2  speed  is  restricted  by  output  time 
which  results  in  the  computation  times  reflected  in  Table  1.  Compu¬ 
tation  times  including  terminal  output  and  error  table  output  only  are 
also  included  in  Table  1.  Clearly,  detailed  flight  history  should  be 
printed  out  only  when  absolutely  necessary. 

Very  high  drag  bodies  such  as  parachutes  cannot  be  treated  by 
this  program. 
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Comparison  of  Computation  Time 
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AVERAGE  TIME  INCREMENT  =  0.  33  SECONDS 
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APPENDIX  I 


E0RTRAN  LISTING 


THE  TEA-2  S0URCE  PROGRAM  IS  WRITTEN  IN  F0RTRAN  LANGUAGE. 
IT  REQUIRES  SINE,  C0SINE,  SQUARE  R00T,  AND  EXP0NENTIAL  LIBRARY 
SUBR0UTINES.  ALS0  REQUIRED  ARE  THE  ATH0SPHERE  (ATM0S)  AND 
table  SEARCH  (TABLE)  SUB-PR0GRAHS  WHICH  ARE  LISTED  AFTER  THE 
MAIN  PR0GRAM.  A  F0RTRAN  LISTING  0F  THE  MAIN  PR0GRAM  F0LL0WS. 


TEA-2  TW0-DIMENSI0NAL,  TW0-STAGE  R0CKET  0R  BALLISTIC  TRAJECT0RY 

(WITH  1959  ARDC  STANDARD  ATM0SPHERE) 

C0MPUTING  and  ANALYSIS  SECTI0N 

TRAJECT0RY  ERR0R  ANALYSIS  ••  TW0  DEGREES  0F  FREED0M 
ANALYSIS  ••  LT.  H.J.  K0PP,PH.D  *  EXTENSI0N  7226A 
PR0GRAMMIN6  ••  MR.  J.N.  NIELSEN  •  EXTENSI0N  73230 

THIS  PR0GRAM  SELECTS  THE  INTEGRATI0N  STEP  LENGTH  SUCH  THAT  THE 
ERR0R  IN  THE  TBTAL  RANGE  IS  LESS  THAN  A  PRESCRIBED  AM0UNt  ANO/0R 
THE  ERR0R  IN  ALTITUDE  IS  LESS  THAN  A  PRESCRIBED  AM0UNT. 


OIMENSI0N  TABL1(2A,2) ,TABL2(2A,2),TABL3(2A,2) ,TRST1(2A,2) , 

1  TRST2(24,2),TMAS(2A,2),TMAS2(24,2) ,NEXT(5) ,TIM(3) , 

2  DWRT(A) ,XD0T(3) ,YD0T(3) ,VD0T(3) ,THO0T(3) 

OIMENSI0N  TN(200) •XN(200) ,YN(2C0) ,THN(200),VN(200),EX(200), 

1  EY(200) .ET(200) ,EV(200) 

C0MM0N  RH0,PRSB,GRAVyCMACH,KDC0N,NATM0S 

F0RMAT  STATEMENTS 


100  F0RMAT  (9I1,3I2,F5.0) 

101  F0RMAT  (8F9.0) 
lOA  F0RMAT  (7F10.0) 

105  F0RMAT  (8F5.0) 

106  F0RMAT  (IHI,A3X,31H  C0MPUTING 


AND  ANALYSIS  SECT I0N/27X ,66H  TW0  DIM 
1ENSI0NAL,  TW0-STA6E  R0CKET  0R  BALLISTIC  TRAJECTORY  PR0GRAM/A1X,37H 
2  *WITH  1959  ARDC  STANDARD  ATMOSPHERE*//) 

(8H1  TIME,7X,6H  RANGE, 6X,9H  ALT  I TUOE, AX, 9H  VELOCITY, AX, 

6H  THETA, 5X,5H  MACH,5X,3H  K0,7X,'»H  THRUST, 6X,7H  WEIGHT, SX, 
DENSITY) 

TIME,7X,6H  RANGE, 6X,9H  ALT  I  TIDE , AX , 9H  VELOCITY, AX, 
THETA, 5X,5H  MACH,5X,3H  K0,7X,7H  THRUST, 6X,7H  WEIGHT, 5X, 
DENSITY) 

(IAX,8H  X-0ERIV,5X,8H  Y-CER I V , 5X , 8H  V-0ERIV,3X,9H  TH-DERIV, 
38X,5H  0RAG,6X,9H  PRESSURE//) 

(F9.3,2FIA.3,F12.3,F11 .3,FIC.A,F9.A,F1A.3,F11.3,E15.5) 
(9X,2FIA.3,F12.3,FII.?,33X,F11.3,E15.5) 

(A6X,20H  END  0F  TR A JECT0RY , I  3  ) 


8H 

(8H 

6H 

8H 


108  FORMAT 
1 
2 

109  FORMAT 
1 
2 

110  FORMAT 
1 

111  FORMAT 

112  format 

113  FORMAT 

115  FORMAT 

llA  FORMAT 

116  FORMAT 

117  FORMAT 

118  FORMAT 
IFF  LAUNCHER) 

121  FORMAT  (A7X,26H 

122  FORMAT  (A5X,30H 

123  FORMAT  (A6X,29H 


(IHL,A6X,26H  COASTING  OF 
(38X,AAH  ACCELERATION  OF 
(1H1,A4X,30H  ACCELERATION 
(1HI,A5X,29H  FREE-FLIGHT 
(1H1,36X,59H  ACCELERATfON 


PAIN  STAGE) 

BOOSTER  AND  MAIN  STAGE) 

OF  MAIN  STAGE) 

OF  MAIN  STAGE) 

OF  BOOSTER  AND  MAIN  STAGE  0 


COASTING  OF 
ACCELERATION 
FREE-FLIGHT 


MAIN  STAGE) 

OF  MAIN  STAGE) 
OF  MAIN  STAGE) 
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126  F0RMAT  (lHly56Xt6H  TEA-2/ 3AX » 52H  TRAJECT0RY  ERR0R  ANALYSIS  **  TUB 
IDEGREES  0F  FREEO0M//) 

127  F0RHAT  (5Xf7lH  ANALYSIS  **  DR.  H.J.  K0PP  *  eXTENSI0N  7226A  •  REP 

10RT  T0  BE  PUBLISHiED/5X,50H  PR0GRAMMING  **  MR.  J.N.  NIELSEN  #  EXTEN 
2SI0N  73230//) 

128  F0RMAT  (8X,103H  THIS  PR0GRAM  SELECTS  THE  INTEGRATI0N  STEP  LENGTH  S 
lUCH  THAT  THE  ERR0R  IN  THE  T0TAL  RANGE  IS  LESS  THAN  A/5X,111H  PRESC 
2RIBE0  AMOUNT  ANO/0R  THE  ERROR  IN  ALTITUDE  IS  LESS  THAN  A  PRESCRIBE 
3D  AMOUNT.  THE  ERR0R  CRITERION  USED  F0R/5X,23H  THIS  TRAJECTORY  WAS 

A  •• ) 

129  F0RMAT  (8X,85H  THE  ABOVE  CONDITION  WAS  SATISFIED,  AND  THE  COMPUTED 
I  ERROR  IN  TOTAL  RANGE  IS  EQUAL  T0,E15.8,9H  PERCENT.//) 

130  FORMAT  (8X,93H  THE  ABOVE  CONDITION  COULD  NOT  BE  SATISFIED,  AS  THE 
ICOMPUTFD  ERROR  IN  TOTAL  RANGE  IS  EQUAL  Ta/6X, E 1 5. 8 , 9H  PERCENT.//) 

131  FORMAT  (28X,39H  THE  ERROR  IN  TOTAL  RANGE  IS  LESS  THAN  ,F6.2, 

1  9H  PERCENT.//) 

132  FORMAT  (2ax,4'iH  THE  ERROR  IN  TERMINAL  ALTITUDE  IS  LESS  THAN  ,F6.2, 

1  29H  PERCENT  OF  MAXIMUM  ALTITUDE.//) 

133  FORMAT  {50X,22H  FINAL  ERROR  EST  IMA TES/57X , 6H  ***♦*) 

134  FORMAT  (5X,7H  RANGE  ,E15.8,14H  RANGE  ERROR  ,E15.8,11H  ALTITUDE  , 

1  E15.8,17H  ALTITUDE  ERROR  ,E15.8/5X,7H  THETA  ,E15.8, 

2  14H  THETA  ERROR  ,E15.8,11H  VELOCITY  ,E15.8, 

3  17H  VELOCITY  ERROR  ,E15.8//) 

135  FORMAT  (51X,19H  INITIAL  C0ND I T  I ONS/ 57X , 6H  *•**•) 

136  FORMAT  (8X,7H  RANGE  , E 15. 8, 3X, ICH  ALTITUDE  ,E15.8,3X,7H  THETA  , 

1  E15.8,3X.10H  VELOCITY  ,E15.8//) 

137  FORMAT  (47X,25H  TABLE  OF  ERROR  ESTIMATES/56X,6H  #*#**) 

138  FORMAT  (1H1,40X,37H  TABLE  OF  ERROR  ESTIMATES  I  CONTINUED ) /56X , 

I  6H  *****) 

139  FORMAT  (11X,5H  TIME,17X,6H  RANGE, 16X,9H  AL T ITUDE , 15X , 6H  THETA, 16X, 

1  9H  VEL0CITY/3OX,12H  RANGE  ERROR , lOX , 15H  ALTITUDE  ERROR, 9X, 

2  12H  THETA  ERROR , lOX , 15H  VELOCITY  ERROR//) 

140  FORMAT  ( 6X , E 15 . 8,4 ( 8X , E 15.8 ) /2 I X, 4 ( 8X , E 15 . 8 ) // ) 

141  FORMAT  (5X,24H  IMPACT  BETWEEN  RANGE  0F,E15.8,9H  FEET  AND,E15.8, 

1  5H  FEET//5X,36H  IMPACT  CENTER  PREDICTED  AT  RANGE  0F,E15.8, 

2  21H  FEET  AND  ALTITUDE  0F,E15.8,5H  FEET) 

142  FORMAT  (8X,90H  THE  ABOVE  CONDITION  COULD  NOT  BE  SATISFIED,  AS  THE 
1C0MPUTEO  ERROR  IN  ALTITUDE  IS  EQUAL  T0/6X , E 15. 8 , 29H  PERCENT  OF  MAX 
2IMUM  ALTITUDE.//) 

143  FORMAT  (8X,91H  THE  ABOVE  CONDITIONS  COULD  NOT  BE  SATISFIED.  THE  CO 
IMPUTED  ERROR  IN  TOTAL  RANGE  IS  EQUAL  T0/3X ,E15. 8,56H  PERCENT,  AND 
2THE  COMPUTED  ERROR  IN  ALTITUDE  IS  EQUAL  T0,E15.8,29H  PERCENT  OF  MA 
3XIMUM  ALTITUDE.//) 

144  FORMAT  (8X,82H  THE  ABOVE  CONDITION  WAS  SATISFIED,  AND  THE  COMPUTED 

1  ERROR  IN  ALTITUDE  IS  EQUAL  T0,EI5.8,11H  PERCENT  0F/5X,17H  MAXIMUM 

2  ALTITUDE//) 

145  FORMAT  (8X,83H  THE  ABOVE  CONDITIONS  WERE  SATISFIED.  THE  COMPUTED  E 
IRROR  IN  TOTAL  RANGE  IS  EQUAL  T0,F15.8,9H  PERCENT , /4X,48H  AND  THE 
2C0MPUTED  ERROR  IN  ALTITUDE  IS  EQUAL  T0,E15.8,29H  PERCENT  OF  MAXIMU 
3M  ALTITUDE.//) 

1  READ  INPUT  TAPE  2,100,  NEXT ( 1 ) , NEXT ( 2) , NE X T ( 3 ) , NEXT (4 ) ,KDC0N, 

1  NTHRST,ICD,M,NXY,  JN,NJR,NS,GAMMA 

GO  TO  (70,74) ,KOC0N 

READ  STATEMENTS 

70  IF  (NEXT(l))  72,72,71 

71  READ  INPUT  TAPE  2,101,  TABLl 

72  READ  INPUT  TAPE  2,101,  TABL2 
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IF  (NEXTO))  74,74,73 

73  READ  INPUT  TAPE  2,101,  TABL3 

74  60  T0  (75,2),NTHRST 

75  IF  (NEXTd))  77,77,404 

404  READ  INPUT  TAPE  2,101,  TRSTl 
60  T0  ( 76, 77) ,M 

76  READ  INPUT  TAPE  2,101,  TMAS 

77  IF  (NEXT! 3))  2,2,78 

78  READ  INPUT  TAPE  2,101,  TRST2 
60  T0  {79, 2), M 

79  READ  INPUT  TAPE 
2  READ  INPUT  TAPE 

IF  (NFXT(l)) 

405  IF  (NEXT(3)) 

406  READ  INPUT 
READ  INPUT  TAPE  2, 

60  T0  408 

407  READ  INPUT  TAPE  2,101,  DWRT ( 4 ) , F ACT ,CKD 

408  WRITE  OUTPUT  TAPE  3,106 


2,101,  TMAS2 

_  2,101,  THDEG0,V0,WCT,DIA,X0,Y0,T£,YF 

405.405.406 

407.407.406 
TAPE  2,101, 


407,406 

, 101,  THRST,D0TM,AREA,BWGT,THR2,O0TN2,AREA2,PRES0 
,101,  TIM,DWRT(1) ,0WRT(2) ,0WRT(4) , CKO, FACT 


CONVERSION  OF  INITIAL  VALLES  TO  CORRECT  UNITS 


T  =  T0 
X  »  X0 

Y  =  Y0 

V  *  V0 

THET  *  THOE60/57. 29577 
THET0  *  THET 
RWCT  =  WGT 

RMASS  =  WCT/32. 17405 
RMA$S0  =  RHASS 
O0TM  *  O0TM/32. 17405 
D0TM2  *  D0TM2/32. 17405 
BMASS  =  BW6T/32. 17405 
CMACH  »  0.0 
OIA  =  OIA/12.0 
DIASQ  >  0IA»0IA 
AREA  »  AREA/ 144. 

PRES0  =  144.O*PRES0 
DWRT(3)  =  OWRT(l) 

IF  (NS)  1001,1001,1002 

1001  NS  *  10 

1002  IF  (GAMMA)  1003,1003,1004 

1003  GAMMA  -  5.0 

1004  IF  (FACT)  1005,1005,1006 

1005  FACT  a  1.0 

1006  60  T0  (5092, 5090 ),NTHRST 

5090  TRSTl (1,1)  *  T0 
TRST1(1,2)  =  THRST 
TRST2(1,1)  =  riM(2) 

TRST2(1,2)  =  THR2 
OELTHR  »  0.1»TIM(l) 

OELTRS  «  0.1*(TIM{3)-TIM(2)  ) 

00  5091  I  *  2,24 

TRSTKI,!)  =  TRSTKI-l,!)  ♦  OELTHR 
TRST1(I,2)  =  THRST 
TRST2(I,1)  »  TRST2(I-l,l)  ♦  OELTRS 
TRST2n,2)  »  THR2 

5091  CRNTINUE 
C 
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C  INITIAL  VALUES  F0R  ERR0R  ANALYSIS  R0UTINE 

C 

5092  AXR  »  1.9E-7 
AYR  »  1.96-7 
ATP  =  2.66-7 
AVTL  »  0.0 
ATTL  *  0.0 
AXTL  *  0.0 
AYTL  =  0.0 
6XLN  =  0.0 
6YLN  *  0.0 
6VLN  =  0.0 
6TLN  =  0.0 
6XNN  »  0.0 
6YNN  »  0.0 
6VNN  =  0.0 
6TNN  =  0.0 
6VN  *  0.0 
6TN  =  0.0 
EXN  *  0.0 
6YN  =  0.0 
Y0U  *  0.0 
0L6Y  =  0.0 
OT0N6  =0.0 
TIME  »  TIM(l) 

NEC  =  0 
NP0S  =  0 
NA  =  0 
NB  =  0 
NTERM  =  1 
NEST  *  I 
NN  »  2 
CNN  *  NN 
N0  =  1 
NSTEP  »  1 
YMAX  *  Y 

IF  (THET)  5991,5991»5992 

5991  JMP  =  2 

60  T0  5993 

5992  JMP  =  1 

5993  JUMP  a  1 

GAM2  =  .01«GANMA 
KPHASE  *  I 
CELT  =  V/32174.05 
OTH  »  OELT 
EPS  =  l.O 
KU  =  I 
KUK  =  1 

ASSIGN  7007  T0  KAP 

ASSIGN  42  T0  NWRT 

G0  T0  (6001f6002»6001) tNXY 

6001  NNXY  «  1 
G0  T0  6003 

6002  NNXY  «  2 

6003  ASSIGN  602  T0  NERR 
C 

60  T0  (80>81)iKD 
60  CF  »  FACT 
G0  T0  82 

81  CF  «  0.39269908*FACT 
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82  J  »  1 
NNN  «  1 
NATP0S  *  I 
88  N  *  1 
LINE  «  2 

ASSIGN  64  T0  NCNT 

INTEGRATieN  C0NSTANtS 

Cl  *  0.62653829  *  CELT 
C2  *  -0.55111241  *  CELT 
C3  *  1.4072559  *  DELT 
C4  =  -0.48268182  *  DELT 
EPSL0N  =  0.001  *  DELT 
WRIT  »  T-EPSL0N 

IF  (NEXT(N) >409, 409*410 

409  NNN  «  2 
G0  T0  50 

410  WRITE  OUTPUT  TAPE  3*114 

61  WRITE  OUTPUT  TAPE  3*109 
WRITE  OUTPUT  TAPE  3*110 

3  DO  45  K  s  1*3 

GO  T0  NCNT* (62*64) 

62  WRITE  OUTPUT  TAPE  3*108 
WRITE  OUTPUT  TAPE  3*110 
ASSIGN  64  TO  NCNT 
CONTINUE 

64  CALL  ATH0S(V*V) 

GO  TO  ( 10*29) (KOCON 
10  GO  TO  (17*25*2l*25)*N 

TABLE  LOOK-UPS  FOR  DRAG  COEFFICIENT 

DURING  FIRST  BURNING  STAGE 

17  CALL  TABLE(TABLl*CMACH*CKn*24) 

GO  TO  29 

DURING  SECOND  BURNING  STAGE 

21  CALL  TABLE(TABL3*CMACH*CK0*24) 

GO  TO  29 

DURING  FREE-FLIGHT 

25  CALL  TABLE(TABL2*CMACH,CKD*24) 

CONTINUE 

29  GO  TO  (95*96) tNTHRST 
95  GO  TO  (200*96*300*96) *N 

THRUST  TABLE  FOR  FIRST  BURNING  STAGE 

200  CALL  TABLE(TRSTl*T*THRST,24) 

GO  TO  (18*96*96)»K 

18  GO  TO  (19*96)*M 

19  CALL  TABLe(TMAS*T*RHASS»24) 

RNASS  «  RNASS/32. 17405 
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60  re  96 


THRUST  TABLE  F0R  SECeND  BURNING  STAGE 

300  CALL  TABLE(TRST2,T,THRST,2A) 

G0  T0  (22,96«96)tK 

22  G0  T0  (23.96) «M 

23  CALL  TABLE(THAS2*T*RHASSt2A) 

RMASS  >  RHASS/32.L7405 
C0NTINUE 

DERIVATIVES 

96  CCD  »  CF»CKD 

DRAG  =  CCD«RH0«V«V«OIASQ 
DOM  «  DRAG/RMASS 
CC0S  =  C0SF(THET) 

CSIN  *  SINF(THET) 

XD0T(K)  *  V*CCl3S 
YD0T(K)  =  V*CSIN 
G0  T0  (501.502.501.502).N 

501  PRESS  =  PRSB*ARGU 

THRUST  =  THRST  *■  AREA*  ( PRESB-PRESS  ) 

VO0T(K)  =  HRUST/RMASS)  -  6RAV«CSIN  -  O0H 
RWGT  =  32. 17A05*RMASS 
G0  T0  503 

502  VO0T(K)  *  -GRAV*CSIN  -  D0M 

503  THO0T(K)  *  >GRAV*CCeS/V 

3A  G0  T0  (35.43,4A),K 

35  G0  T0  NMRT,(36,39»41.42} 

36  IF(OWRT(N))  37,37,38 

37  ASSIGN  41  T0  NWRT 
G0  T0  41 

38  ASSIGN  39  T0  NMRT 
CONTINUE 

39  IF(WRTT-T)  40,40,42 

40  WRIT  *  T+OWRT(N)-EPSL0N 

41  THDEG  *  57.295?7*THET 
THDDT  =  57.29577*THD0T( 1) 

WRITE  0UTPUT  TAPE  3,111,  T , X , Y , V, THDEG , CMACH , CCD, THRUST ,RHGT,RH0 
WRITE  0UTPUT  TAPE  3,112,  XCBTI I ) ,YC0T( 1 ) ,VC0T ( 1 ) ,THDDT, DRAG, PRESS 
LINE  *  LINE  +  1 
IF  (28-LINE)  63,63,  42 
63  ASSIGN  62  T0  NCNT 
LINE  *  I 
CONTINUE 

M0OIFIEC  THIRD  0RDER  RUNGE-KUTTA  INTEGRATION 

42  X  =  X  ♦  Cl*XD0T(K) 

Y  =  Y+C1*YD0T(K) 

V  *  V+C1*VD0T(K) 

THET  *  THET+C1*THO0T(K) 

GO  TO  45 

43  X  *  X  ♦  C2*XD0T(K) 

Y  *  Y+C2*YD0T(X) 

V  *  V+C2*VD0T(K) 

THET  »  THET4C2*THD0T(K) 


50 


non  o  o  o  o 


G0  re  A5 

AA  X  =  X  +  C3*XC0T(K-1)  ♦  CA#XO0T(K) 

Y  =  Y+  C3»YO0T(K-l  )+CA«YDClT(K) 

V  =  V+  C3*VO0T(K-l )+CA»VC0T(K) 

THET  *  THET+C3»THD0T(K-1)+CA*THC0T(K) 

T  =  T+OELT 

RMASS  =  RMASS-O0TM*CELT 
A5  C0NTINUE 

M.K0PP  ERR0R  ANALYSIS  R0UT INE  ISIHPLIPIEO) 

NSTEP  *  NSTEP  ♦  I 
G0  T0  (599,599,92,92) ,KU 
599  IF  (NN  -  NSTEP)  601,601,92 

601  G0  T0  (1311,92),NTERN 

1311  IF  {Y-Y0U)  92,1312,1312 

1312  G0  T0  NERR, (602,603,605) 

602  ASSIGN  603  T0  NERR 
NSTEP  *  NN  -  I 

VV  =  V 
THTiH  »  THET 
XX  =  X 
YY  =  Y 
RHT  =  RMASS 
TT  »  T 
ODTT  »  DELT 
G0  T0  3 

603  ASSIGN  605  T0  NERR 
ASSIGN  A2  T0  NWRT 
NSTEP  =  NN  -  2 

OH  *  DELT 
VH  »  V 
THH  *  THET 
XH  =  X 
YH  =  Y 

C0MPUTATI0N  0F  PR0PA6ATieN  ERR0R  COEFFICIENTS 

C  =  (2.O*RH0*CCD*OIASQ*V*V)/(RMASS»6RAV) 

CSQ  =  C*C 

PC  =  C0SF(THET) 

PS  =  SINF(THET) 

PS  «  A8SF(PS) 

AI3  *  PC»(1.03  -  .01A*C  ♦  .C003*CSC) 

AlA  =  1.2*PS  ♦  .28  -  .OOOUC 

A23  *  (PS<-.O8)*(l.0-.02*C*.C00A«CSC) 

AAA  *  1.2*PS  ♦  .21  ♦  .0003«C 
PCP7  *  PC  ♦  .07 
A2A  «  AAA«PCP7 

A33  *  .019  ♦C*(“l.l  ♦  C*(.0A  -  .0002*0' 

A3A  s  (PC+.l)*(  1.02+. 018*0 
AA3  *  PCP7*(1.A-.02*C+.0003*CSQ) 

DELTA  =  A33*AAA  -  A3A*AA3 

AVR  *  (THRUST/RW6T)*2.9E-7  ♦  r*3.3E-7  ♦  2.0E-7 
AFT  »  A33-AAA 

ARGM  «  AFT*AFT  ♦  A.0*A3A*AA3 

EPS  «  0.5*(A33  *  AAA  *  SORTF(ARGM)) 
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GKK  «  (EPS-A33)/A34 
EPS  *  GRAV«EPS/V 

V  =  VV 
THET  =  THTH 
X  =  XX 

Y  =  YY 
T  =  TT 
RMASS  *  RMT 

DELT  »  0.5*DELT 

604  Cl  =  0.62653829»DELT 
C2  =  -0.5511l24l*0ELT 
C3  *  1.4072559#0ELT 
C4  =  -0.48268182*DELT 

G0  T0  3 

605  EH  =  EPS*DH 
ASSIGN  602  T0  NERR 
CLAM  =  1.0  ♦  EH 
EH02  *  0.5*EH*EH 

C0MPUTATI0N  0F  TRUNCATI0N  ERR0R  COEFFICIENTS 

A  *  ABSF(VH-V)  -  0H*( ITHRUST/RMASS)*11.0E-8  ♦  GRAV*7.0E-8  ♦  D0M« 
1  19.0E-8) 

IF  (A)  5001,5002,5002 

5001  A  »  0.0 

5002  8  =  A8SF(THH-THET)  -  ( GRAV«OH« 1 .OE-8 ) /V 
IF  (8)  5003,5004,5004 

5003  B  -  0.0 

5004  DH3  =  DH**3 

CC  *  l.75*V*0H3 

AVT  *  (A*ABSF(A33)  ♦  VH»A34«B)/CC 
IF  (AVT-AVTL)  4002,4001,4001 

4001  AVTL  *  0.8*AVT 

G0  T0  4003 

4002  AVT  *  AVTL 
AVTL  *  0.8»AVT 

4003  CONTINUE 

ATT  *  (A43*A  +  VH*A44*B)/CC 
IF  (ATT-ATTL)  4005,4004,4004 

4004  ATTL  *  0.8*ATT 

G0  T0  4006 

4005  ATT  a  ATTL 
ATTL  *  0.8*ATT 

4006  CONTINUE 

PHB  a  A44*AVT-A34*ATT 

PH6B  a  A44*AVR  -  A34*ATR 

PSB  a  -A43*AVT  ♦  A33*ATT 

PSBB  a  -A43*AVR  ♦  A33«ATR 

EV0V  a  EVN/VH 

VH4  a  VH*DH*0H3 

TH3  a  2.0*DH3 

QQ  a  1. 142857 1/VH4 

DHV7  a  0H*V»7.0E-8 

DX2  a  ABSF(XH-X)  -  0HV7 

IF  (0X2)  5005,5006,5006 

5005  0X2  a  0.0 

5006  AXT  a  QQ«0X2 

IF  (AXT-AXTL)  4008,4007,4007 

4007  AXTL  a  0.8»AXr 
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non 


G0  T0  4009 

4008  AXT  =  AXTL 
AXTL  »  0.8»AXT 

4009  CONTINUE 

0V2  =  ABSF(YH-Y)  -  DHV7 
IF  (DY2)  5007,5008,5008 

5007  0Y2  *  0.0 

5008  AYT  =  QQ*0Y2 

IF  (AYT-AYTL)  4011,4010,4010 

4010  AYTL  *  0.8*AYT 
G0  T0  6005 

4011  AYT  =  AYTL 
AYTL  a  0.8*AYT 

TRAJECTORY  ESTIMATES  F0R  THRUST  PHASES 

6005  OYK  =  Y 
QYH  =  Y 
QXK  a  X 
QXH  a  X 
QVK  a  V 
QVH  a  V 
QTHK  a  FHET 
QMH  a  RMASS 
OMD0T  a  DOTH 

60  T0  (6006, 6026, 6008, 6028), NEST 

6006  00  6366  II  a  1,30 

IF  (TRSTKII,!)  -  T)  6366,6166,6266 
6166  QTK  a  TRSTK  II  ,2) 

QTH  a  TRSTK  1 1*1, 2) 

QMH  a  TMAS(  IKl, 21/32. 17405 
C0NT  a  TRSTK  I  Kl,  1)  -  TRSTKII,!) 

60  TO  6466 

6266  OIK  a  TRSTKII-1,2)  ♦  K  T-TRST  K  I  I-l » 1 ) )  /  (  TRSTl  ( 1 1 , 1  )-TRSTl  ( I  I-l  tl 
1)))  ♦  (TRSTK  I  I,2)-TRSTK  I  1-1,2)) 

OTH  a  TRSTKII, 2) 

QMH  a  TMAS( 11,21/32.17405 
C0NT  a  TRSTK  11,1)  -  T 
60  T0  6466 
6366  CONTINUE 
6466  KLER  »  i 

60  T0  ( '^10,6011  ),M 
6008  DO  6388  II  =  1,30 

IF  (TRST2(II,L)  -  T)  6388,6188,6288 
6188  QTK  a  TRST2( II ,2) 

QTH  a  TRST2( 11+1,2) 

QMH  a  TMAS2( I  1  +  1,21/32.17405 
C0NT  a  TRST2( I  I  +  l,  1)  -  TRST2(II,1) 

60  T0  6488 

6288  QTK  =  TRST2( I I-l ,2  )  ♦  K r-TRST2 ( I I-l ,  1 ) ) / ( TRST2 ( 1 1 , 1 )-TRST2 ( I I-l , 1 
1)))*  (TRST2(II,2)  -  TRST2( 11-1,2) 1 
QTH  a  TRST2( I  I ,2) 

QMH  a  TMAS2I II ,21/32.17405 
C0NT  a  TRST2(II,1)  -  T 
60  TO  6488 
6388  CONTINUE 
6488  KLER  a  2 

OMOOT  a  O0TM2 
60  TO  (6010, 6011), M 
6010  QM  a  (QMH-RMASS)/C0NT 
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G0  T0  6012 

6011  QM  =  -DMD0T 

QMH  =  RMASS  -  DMD0T«C0Nr 

6012  QT  =  (QTH  -  THKST)/C0NT 
ALPHA  =  GRAV»SINF(THET) 

BETA  =  RH0*OI ASQ*QVK*CKD«CF 
AQT0B  =  ALPHA  +  CT/BETA 
QEX  =  BETA/OM 

UN  =  QM  4-  BETA 

QVH  =  QTH/BETA  -  ( CMH*AQT0B/UN )  ♦  (QVK  ♦  AG T0B«RMASS/UK  -  THRST/ 

1  BETA)«{ (RMASS/QMH)««CEX) 

PT3V  =  O.A*(QVH+QVK  )*C0NT 
CSKAY  =  C0SF(QTHiK) 

QXK  =  QXK  ♦  PT3V*CSKAY 
QYK  =  QYK  +  PT3V»SINF{UTHK  ) 

QTHK  *  GTHK  -  2.0»6RAV*CSKAY*C3NT/ ( CVK  +  CVH  ) 

QVK  =  QVH 
II  =  II  +  1 

6013  00  6021  I  =  1 1 f 30 
CALL  ATM0S(QYK,QVK) 

GO  T0  (6014,6016) tKLER 

6014  QTK  =  QTH 

QTH  =  TRSTKI.Z) 

C0NT  »  TRSTKI.l)  -  TRSTl(I-l,l) 

TM0UT  =  TRSTK  1,1) 

60  T0  (6614,6615) ,KOC0N 

6614  CALL  TABLE (TA8Ll,CMACH,CKC, 24) 

6615  G0  T0  (6015, 6019), M 

6015  QMK  =  QMH 

QMH  =  TMAS( I,2)/32.17405 
G0  T0  6018 

*6016  QTK  =  QTH 

QTH  =  TRST2(I,2) 

C0NT  »  TRSr2(I,l)  -  TRST2(I-l,l) 

TM0UT  =  TRST2(  1,1) 

G0  T0  (6616,6617) ,KOCBN 

6616  CALL  TABLE(TA8L3,CMACH,CKD,24) 

6617  60  T0  (6017,6019) ,M 

6017  QMK  =  QMH 

QMH  =  TMAS2( I,2)/32.17405 

6018  QM  =  (OMH-QMK)/C0NT 

G0  T0  6020 

6019  QM  =  -OMC0T 
QMK  =  QMH 

QMH  =  QMK  -  DMr0T»C0NT 

6020  CT  =  (QTH-QTK) /C0NT 
SNKAY  =  S  INF (QTHK) 

CSKAY  *  C0SF(QrHK) 

ALPHA  *  GRAV*SNKAY 

BETA  =  RH0*OI ASQ«QVK*CKC»CF 
AQT0B  =  ALPHA  ♦  QT/BETA 
UN  *  QM  +  BETA 
QEX  =  BETA/QM 

QVH  *  (QTH/BETA)  -  (QMH*ACT08/UN )  ♦  (QVK  +  ( AOTBB*OMK/LN )  -  (QTK/ 
1  BETA) )*( (QMK/CMH)**QEX) 

PT3V  =  O.4»(QVH+QVK)*C0NT 
QXK  =  QXK  ♦  PT3V*CSKAY 
QYK  »  QYK  ♦  PT3V»SNKAV 

QTHK  *  QTHK  -  2.0»GRAV«CSK AY*CBNT/ ( CVK+QVH ) 
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o  o  o  o  o  o 


QVK  =  QVH 

IF  (TIME-TM0UT-C0NT*l.OE-'!»)  6022f  6022t  6021 

6021  CONTINUE 

6022  Gi3  T0  (6024t6028)  »KLER 
6026  IF  (NEXT(2) )6028, 6028, 6025 

TRAJECTORY  ESTIMATE  F0R  INTERMEDIATE  BALLISTIC  PHASE 

6025  OMK  =  QMH  -  BMASS 
C0NT  =  TIM(2)  -  TIM(l) 

G0  T0  6027 

6026  QMK  =  RMASS 

C0NT  =  TIM(2I  -  T 

6027  PT6V  =  O.6*QVK#C0NT 
CSKAY  =  C0SF(grHiK) 

SNKAY  =  SINF(grHiK) 

QXK  =  QXK  ♦  PTftV*CSKAY 

QYK  =  QYK  ♦  PT6V*SNKAY 

QTHK  =  QTHK  -  {GRAV«CSKAY/C VK > *C0NT 

G0  T0  (6627, 6628), KCC0N 

6627  CALL  TABLE ( TABL2,CMACH ,CKC , 26 ) 

6628  QVK  =  QVK  +  (-GRAV»CSKAY  -  RH0*CF*CKD«QVK*CVK*D lASQ/QMK )»CCNT 
KLER  =  2 

II  =  2 

DMO0T  =  D0TM2 
CTH  *  TRST2( 11-1,2) 

C0  T0  6013 

TRAJECTORY  ESTIMATE  FOR  TERMINAL  BALLISTIC  PHASE 

6028  U  =  SINF(GTHK) 

QD0TY  *  QVK*U 

TMTN  =  (OO0TY  ♦  SQRTF( QDBTY«QD0TY  ♦  2. 0*GRAV»CYK ) ) /GRAY 
G0  T0  (6129,6128,6129) ,NXY 

6128  G0  T0  (6129,5016) ,JMP 

6129  CSK  =  C0SF(QTHK) 

QOOTX  =  0.6*GVK*CSK 
PLUS  -  1.0/(QVK*CSK) 

G0  T0  (6227, 6228), KOC0N 

6227  CALL  ATM0S ( QYK ,QVK ) 

CALL  TABLE(TA0L2,CMACH,CKC,26) 

6228  BET  =  0 .0023769*0 1 ASQ»eVK#CKD»CF/ 32 . 1 7605 
SN0CS  =  U/CSK 

DU  =  0.05*(U+0.9995) 

DU0G  =  OU/32. 17605 
C 

6229  IF  (QYK-5.8E6)  6029,6030,6031 

6029  B0G  =  BET 

G0  T0  6063 

6030  B0G  =  0.1*BET 

G0  T0  6063 

6031  IF  (QYK-1.0E5)  6030,6032,6033 

6032  B0G  =  0.01*BET 

G0  T0  6063 

6033  IF  (QYK-1,6E5)  6032,6036,6035 
6036  B0G  =  0.00l*BET 

G0  T0  6063 

6035  IF  (QYK-2.3E5)  6036,6036,6037 

6036  B0G  «  BET*1.0E-6 

G0  T0  6063 
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o  n  o 


<037  IF  (QYK-3.0F5)  6036,6038t 6C39 
6038  B3C  =  BET*l.0E-6 
00  T0  60A3 

<039  IF  {QYK-4.0F5)  60 38 ,60A0, 6C^ I 

6040  B0G  s  BFT*l.0E-8 

GB  T0  6043 

6041  IF  (QYK-6.0E5)  6040 , 6042 1 6042 

6042  B3G  =  BET*1.0E-13 
C 

6043  AA  =  B0G*SN0CS  ♦  PLUS 
0MU  =  (  1.0+U)*(l.0-U) 

R30T  =  SQRTF(0MtJ» 

AU  =  (R00T*AA  -  B3G*U)**2 
C 

OO0TY  =  U/(2MU*AU) 

OYH  =  QYK  ♦  O.H*QD0TY*DU0G 
OO0TX  =  1.0/ ( AU*R00T) 

QXH  =  QXK  ♦  QDOTX*nU0G 
G0  T0  (6444,6046) ,JPP 
6444  G0  T8  ( 6044 , 6046 ), JUMP 

6044  IF  (QYH  -  QYK)  6045,6045,6046 

6045  YMAX  *  QYK 
JUMP  *  2 

6046  IF  (QYH-YF)  6049,6048,6047 

6047  QYK  =  QYH 
QXK  -  QXH 
U  =  U  -  DU 

IF  (l.ChFU)  6048,6048,6229 
C 

6048  XMAX  s  QXH 

G0  T0  5014 

6049  XMAX  *  QXK  ♦  ( QXH-QXK ) • ( QYK-YF  ) / ( OYK-QYH ) 

ESriMAFES  0F  ERRORS  0VER  N  STEPS 
5014  JUMP  *  1 

IF  (NN-NSTEP)  8002.8002,8001 

8001  CNN  =  NSTEP 

8002  CLAMN  =  CLAM**CNN 

GL  =  (CLAMN-l.0)/{CLAM-1.0) 

VHH  =  VH»DH 

PHH  =  ABSF(PH8«0H3  ♦  PH08) 

PSH  =  ABSF(PS0»OH3  ♦  PSBB) 

IF  (PHH-PSH)  4013,4014,4014 

4013  CONE  *  PSH/ABSF(OELTA) 

G0  T0  4015 

4014  C0NE  =  PHH/ABSF( DELTA) 

4015  ET0K  =  ETN/ABSF(6KK) 

IF  (EV0V-ET0K)  4016,4017,4017 

4016  CTW0  *  ET0K 
G0  T0  4018 

4017  CTW0  =  EV0V 

4018  CTHR  *  A13  *  6KK*A14 

CF0R  «  A23  *  GKK«A24 

C12  »  C0NE  +  CTW0 

EVLN  »  V*I-C0NE  ♦  CLAMN*C12) 

ETLN  =  GKK*EVLN/V 

EXLN  *  EXN  *  VHH*(CNN*(AXT«0H3  *  AXR)  -  C0NE*CTHR  *  CTHR«C12*GL) 

EYLN  >  EYN  *  VHH*(CNN*( AYT«0H3  ♦  AYR)  -  C0NE*CF8R  *  CFeR*C12*6L) 


o  o  o 


8003  EVN  < 
ETN  * 
EXN  » 
EYN  * 
G0  T8 


IF  (NN-NSTEP)  8003»8003,80C4 


0.5*IFVLN  ♦  EVNN) 
0.5*(ETLN  ♦  ETNN) 
0.5«(EXLN  ♦  EXNN) 
0.5*(EYLN  +  EYNN) 
8005 


8004  EVN  « 
FTN  « 
EXN  = 
EYN  s 

8005  G0  T0 


EVLN 

ETLN 

EXLN 

EYLN 

(7111,7222.7333) fKPHASE 


FRR0R  ANALYSIS  INITIALIZATION  F0R  START  0F  PHASE 

7222  KPHASE  =  3 
NN  =  4 
NSTEP  =  2 
CNN  *  2.0 
AVTL  =  0.0 
ATTL  =  0.0 
AXTL  =  0.0 
AYTL  =  0.0 
EVNN  *  0.0 
ETNN  »  0.0 
EXNN  =  0.0 
EYNN  =  0.0 
EXAVE  =  EXN 
EYAVE  =  EYN 
EVAVE  =  EVN 
ETAVE  »  ETN 
T00  =  T 
X00  *  X 
Y00  *  V 
V00  =  V 
TH0C  =  THET 
RM00  =  RMASS 
OTH  »  V/32174.05 
ASSIGN  42  T0  NWRT 
G0  T0  637 
C 

7333  KPHASE  =  I 
T  =  T00 
X  =  X00 

Y  =  Y00 

V  =  V00 
THET  =  TH00 
RMASS  =  RM00 
EXN  =  EXAVE 
EYN  =  EYAVE 
EVN  *  EVAVE 
ETN  =  ETAVE 
ASSIGN  36  T0  NHRT 
G0  T8  609 

C 

Till  IF  (NN-2)  7001,7001,7002 
7001  T  =  T0 
X  =  X0 

Y  =  Y0 

V  *  V0 
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non 


THET  ^ 

=  THETO 

RMASS 

s  RMASS0 

EXN  = 

0.0 

EYN  = 

0.0 

EVN  = 

0.0 

ETN  = 

0.0 

7002  CNN  = 

NS 

NN  =  NS 
NSTEP  =  I 
ASSIGN  36  Td  NMRT 
IF  (SENSE  SWITCH6)  300A,30C6 
300A  CALL  PDUMPIZF, A13,  1) 

3006  CONTINUE 

ESTABLISHING  STARTING  VALUES  FOR  DOUBLE  FALSE  POSITION 

609  KLK  =  1 

DT  =  0.035*VH/GRAV 

610  CLAM  =  1.0  ♦  fcPS*DT 
GNN  =  TMTN/OT 
CLAMN  =  CLAM*»GNN 

GL  =  (CLAMN-1.0)/(CLAM-1.0) 

0T3  =  0T**3 

PHT  =  A8SF(PHa»DT3  ♦  PHBB ) 

PST  =  ABSF|PSB*DT3  ♦  PSBB) 

IF  (PHT-PST)  4020,A021,AC2l 
A020  CONE  =  PST/A8SF(0ELTA) 

GO  TO  4022 

4021  CONE  =  PHT/ABSF(CELTA) 

4022  ET0K  =  ETN/ABSF(GKK) 

EV0V  =  EVN/VH 

IF  (EVOV-ETOK)  4023f4024,4C24 

4023  CTWO  =  ET0K 
GO  TO  4025 

4024  CTWO  *  EV0V 

4025  C12  =  CONE  ♦  CTWO 
VHH  =  VH»CT 

GO  TO  (4026,4027) ,NNXY 

4026  EXNN  *  EXN  +  VHH* ( GNN* ( AXT •CT3  ♦  AXR )  -  CONE*CTHR  ♦  CTHR*C12*GL) 
XIOT  =  EXNN  -  GAM2*XMAX 

GO  TO  4028 

4027  EYNN  =  EYN  ♦  VHH* ( GNN* ( AYT*CT3  ♦  AYR)  -  CONE*CF0R  ♦  CFBR*C12*GL) 
XI DT  =  EYNN  -  GAM2*YMAX 

C 

4028  GO  TO  (611 ,615,615,615,617) , KLK 

611  IF  (XIOT)  612,613,613 

612  DTH  =  DT 

GO  TO  630 

613  Xin2  =  XIDT 
DT2  =  DT 

614  DT  *  0.1*DT 
KLK  =  KLK  ♦  1 

GO  TO  610 

615  IF  (XIDT)  616,614,614 

616  XIDl  *  XIDT 
DTI  *  DT 

GO  TO  618 

617  GAM2  »  2.0*GAM2 

GO  TO  609 
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c 

c 

c 

618 


4030 

40  31 

4032 

4033 

4034 

4035 

5020 


5321 


5022 


660 

5025 

5026 

5027 


5028 

5029 


661 

619 


620 


621 

622 


623 


C 


DOUBLfc  FALSE  P0SITI0N  R0UT1NE  F0R  F(DELT) 

PT0  =  (0T1*XI0?  -  DT2*XIDn/{X  102  -  XIOl) 

GNN  =  TMTN/DT0 
CLAM  =  1.0  +  EPS»DT0 
CLAMN  =  CLAM**GNN 
DT3  =  DT0»*3 
VHH  *  VH*DT0 

GL  =  (CLAMN-1.0)/(CLAM-1.0) 

PHT  =  ABSF(PHd*DT3  ♦  PHB8) 

PST  =  ABSF(PSB*0T3  ♦  PSBB) 

IF  (PHT-PST)  4030,4031,4031 
C0NE  =  PST/ABSF(CELTA) 

G0  T0  4032 

C0NE  =  PHT/A8SF(DELTA) 

ET0K  =  ETN/ABSF(GKK» 

IF  (EV0V-ET0KI  4033,4034,4034 

CTW0  =  ET0K 

G0  T0  4035 

CTW0  =  EV0V 

C12  =  C0NE  +  crw0 

G0  T0  (5020,5021) ,NNXY 

EXNN  =  EXN  +  VHH*(GKN*{ AXT*CT3  ♦  AXrt)  -  C3NE*CTHR  +  CTHR*C12»GL) 
Xin0  =  EXNN  -  GAM2*XMAX 
G0  T0  5022 

EYNN  =  EYN  ♦  VHH* ( GNN* ( AYT«CT 3  ♦  AYR)  -  C0NE*CF0R  +  CFER*C12*G\) 
Xin0  =  EYNN  -  GAM2*YMAX 
DIFH  =  (0T1-CT2)/DT2 
IF  (0.05-ABSF(QIFH) )  661,660,660 
G0  r0  (5025, 5026), NNXY 
G0  70  (5029,5029,5027) ,NXY 
60  T0  (5029,5029,5028) ,NXY 
KLK  a  1 
OT  =  nT0 
NNXY  *  2 
G0  T0  610 
NNXY  =  1 
DTH  s  OT0 
NEG  =  0 
NPeS  *  0 
G0  T0  630 

IF  (XIO0)  619,620,620 
NEG  *  NEG  ♦  1 
NP0S  *  0 

60  70  (621,626) , NEG 
NPBS  =  NP0S  ♦  1 
NEG  =  0 

G0  70  (621, 627), NPES 
IF  (XIO0)  622,623,623 
on  *  070 
XIOl  =  XIO0 
G0  70  618 
072  *  070 
XID2  *  XIO0 
G0  70  618 


626  IF  (XIOl)  627,628,628 

627  XI02  «  XI02/2.0 
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NEG  *  0 
NPeS  *  0 
G0  re  618 

628  XlDl  »  XIOl/2.0 
NEG  *  0 

NPBS  *  0 
G0  T0  618 

629  IF  (XIDl)  628,627i627 
C 

630  G0  T0  (632f632t632»631 I tN 

631  NN  »  NS 
CNN  a  NS 
KPHASE  =  1 
NSTEP  =  1 

ae  T0  636 

632  CNN  a  NS 

IF  (CNN*CTH+T  -  TIM(N))  636,636t633 

633  CNN  =  ( riM(N)-T)/DTH 
NN  a  CNN 

CCN  a  NN 

IF  (CNN-CCN)  63A»634t63S 
636  CNN  a  CCN 

G0  T0  6355 

635  CNN  a  CCN  ♦  1.0 

6355  NSTEP  =  1 
NN  a  CNN 

6356  NEST  *  NEST  ♦  1 

IF  (NEXT(NEST))  6356, 6356 »6 357 

6357  KPHASE  *  2 

OTH  a  (TIM(N)-T)/CNN 

636  CLAMN  a  CLAM»»CNN 

GL  a  (CLAKN-l.OI/ICLAH-l.O) 

EVNN  *VH*(-C0NE  ♦  CLAHN*C12) 

ETNN  a  GKK»EVNN/VM 

EXNN  a  exN  ♦  VHH»(CNN»(AXT«CT3  ♦  AXR)  -  r.0NE*CTHR  ♦  CTHR*C12»GL) 
EYNN  a  eYN  ♦  VHH*(CNN*(AXT*DT3  ♦  AYR)  -  C0NE*CF0R  ♦  CF0R«C12*GL) 
C 

G0  T0  (7006,7003,7006,637) ,KU 
7003  Y0U  a  yF  +  EYN 

OLEY  a  (EYNN  -  EVN)/CNN 
EHRA  a  -OLEY 
7006  EV(N0)  a  EVN 

ET(N0)  a  ETN*57. 29577 

EX(N0)  a  EXN 

EY(N0)  a  EYN 

TN(N0)  a  tT  ♦  OH 

XN(N0)  a  XH 

YN(N0)  a  YH 

VN(N0)  a  VH 

THN(N0)  a  THH  *  57.29577 
N0  a  N0  >  1 

637  DELT  =  OTH 

Cl  a  0.62653829»0ELT 
C2  a  -0.55111261*0ELT 
C3  a  l,4072559*DELT 
CA  a  -0.4B268182»0ELT 
G0  T0  (92, 92, 69), KPHASE 
C 

69  THOEG  a  57.295 r7«THET 

THODT  a  57.29577*THO0T(3) 

60 


o  o  n 


RWGT  *  32.17405  *  RMASS 

WRITE  0UTPUT  TAPE  3, 111,  T  ,  X ,Y , V, THOEG.CMACH , CCD, THRUST  , RWGT, RH0 
WRITE  0UTPUT  TAPE  3,112,  XC0T( 1 ), YC0T( 1 ) ,VC0T ( 1 ), THOOT , DRAG, PRESS 

50  N  =  N+1 

IFINEXTINJ)  50,50,51 

51  G0  T0  (3, 55, 56, 59, 7020), N 

55  THRUST  *  0.0 
D0TM  =  0.0 
AREA  =  0.0 

RMASS  =  RMASS-OMASS 
RM00  =  RMASS 
RWGT  »  32.17405  *  RMASS 
G0  T0  (801,802),NNN 

801  WRITE  0UTPUT  TAPE  3,115 
G0  T0  999 

802  WRITE  0UTPUT  TAPE  3,121 
NNN  =  I 

G0  T0  999 

56  THRST  =  THR2 
n0TM  =  O0TM2 

AREA  =  AREA2/144.0 
KASE  =  2 

G0  T0  (803,804) , NNN 
903  WRITE  OUTPUT  TAPE  3,116 
G0  T0  999 

804  WRITE  OUTPUT  TAPE  3,122 
NNN  >  1 

G0  T0  999 

59  thrust  =0.0 
nOTM  =  0.0 
AREA  =  0.0 
GO  T0  (805,806), NNN 

805  WRITE  OUTPUT  TAPE  3,117 
GO  TO  999 

806  WRITE  OUTPUT  TAPE  3,123 
NNN  *  1 

999  EPSL0N  *  O.OOl»OELT 
NEST  =  N 
TIME  =  TIM(N) 

ASSIGN  42  TO  NWRT 
ASSIGN  64  TO  NCNT 
LINE  =  2 

WRITE  OUTPUT  TAPE  3,109 
WRITE  OUTPUT  TAPE  3,110 
GO  TO  3 

TEST  FOR  END  OF  TRAJECTORY 

92  GO  TO  KAP, ( 7007,7010,7013) 

7007  IF  (YOOTIU)  7008,3,3 
1008  KU  *  2 

ASSIGN  7010  TO  KAP 
YMAX  =  Y 
JMP  =  2 

7010  IF  (Y-Y0U)  7012,7011,7011 

7011  YOU  =  YOU  ♦  DLEY 


j1 


YLST  »  Y 
60  T0  3 

7012  ASSIGN  42  T0  NWRT 
NTERM  *  2 

Y0NE  »  Y 
YTW0  *  YLST 
DTTW0  =  -CELT 

DgLT3  =  -OTTW0«{Y0NE-Y0U)/(YTW0-Y0NE) 

DELT  =  DELT3 
ASSIGN  7013  T0  KAP 
60  T0  604 

7013  IE(Y-Y0U)  7014,7015,7015 

7014  Y0NE  *  Y 
.N3  =  0 

NA  *  NA  ♦  1 

IF  (NA-2)  6061,6080,6080 

6080  YTW0  »  Y0U  ♦  O.5*(YTW0-Y0U) 

6081  OT0NE  =  0ELT3 

0ELT3  »  ( ( Yrw0-Y0U)*DT0NE  -  I Y0NE-Y0U ) *0TT W0 ) / I YTW0^Y0NE ) 
DELT  =  0ELT3  -  DT0NE 
60  T0  7016 

7015  YTW0  =  Y 
NA  =  0 

NB  =  NB  ♦  I 

IF  (NB-2)6083, 6082, 6082 

6082  Y0NE  =  Y0U  +  0 .5*  ( YBNE-Y0ij  ) 

6083  OTTVj0  *  0ELT3 

DELT3  *  {  (  YTW0-YBli)»DT0NE  -  I  Y0NE-Y0U )  *DTT W0 ) /  I YTW0-Y0NE ) 
DELT  =  0ELT3  -  OTTW0 

7016  IF  (YTW0-YZNE-.Ol»EYN)  7017,7017,604 

7017  60  T0  (7018,7018,7019,7020) ,KU 

7018  XUN0  =  X 
KU  *  3 
NN  *  NS 
CNN  *  NSTEP 

nstep  =  1 

TT  *  T 
XH  =  X 
THH  =  THET 
YH  =  Y0U 
Y0U  *  YF 
DLEY  *  0,0 
VH  =  V 
YLST  s  Y 
OT0NE  =  0.0 
DH  =  0,0 

ASSIGN  7010  T0  KAP 
60  T0  6005 

7019  EV(N0)  *  EVN 

ET(N0)  *  ETN»57. 29577 

EX(Ne)  =  EXN 

EY(N0)  *  EYN 

TN(N0)  =  T 

XN(N0)  »  X 

YNINE)  *  Y 

VN(N0)  -  V 

THN(N0)  *  THET  *  57.29577 
KU  a  4 
NN  a  NS 
CNN  *  NSTEP 
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NSTEP  *  I 

Tr  =  T 

XH  a  X 

THH  *  THET 

YH  =  Y0U 

VH  =  V 

YLST  *  Y 

DH  =0.0 

Y0U  =  YF-gYN 

OLEY  s  ER«A 

DT0NE  =  0.0 

ASSIGN  7010  T0  KAP 

THDEG  =  THET*57. 29577 

WRITE  0UTPUT  TAPE  T , X, Y , Vf THDEG .CMACH ,CCO , THRUST ♦RWGT , RH0 

401  WRITE  0UTPUT  TAPE  3,113,  NJR 
WRIT  =  T  +  1000. 

ASSIGN  42  T0  NWRT 
G0  T0  6005 

7020  LINE  =  27 
XO0S  =  X 

WRITE  OUTPUT  TAPE  3,126 
WRITE  OUTPUT  TAPE  3,127 
WRITE  OUTPUT  TAPE  3,128 
GO  TO  (5051,5053,5052) ,NXY 

5051  WRITE  OUTPUT  TAPE  3,131,  GAKMA 
PERX  =  1OO.*EX(N0)/XN(N0) 

GO  TO  5054 

5052  WRITE  OUTPUT  TAPE  3,131,  GARMA 

5053  WRITE  OUTPUT  TAPE  3,132,  GAKMA 
PERY  =  1OO.»EY(N0)/YMAX 

5054  IE  (GAMMA-100. *GAM2)  5055,5056,5056 

5055  GO  TO  (5551,5552,5553) ,NXY 

5551  WRITE  OUTPUT  TAPE  3,130,  PERX 
GO  TO  5057 

5552  WRITE  OUTPUT  TAPE  3,142,  PERY 
GO  TO  5057 

5553  WRITE  OUTPUT  TAPE  3,143,  PERX, PERY 
GO  TO  5057 

5056  GO  TO  ( 5554,5555,5556) ,NXY 

5554  WRITE  OUTPUT  TAPE  3,129,  PERX 
GO  TO  5057 

5555  WRITE  OUTPUT  TAPE  3,144,  PERY 
GO  TO  5057 

5556  WRITE  OUTPUT  TAPE  3,145,  PERX, PERY 

5057  WRITE  OUTPUT  TAPE  3,133 

WRITE  OUTPUT  TAPE  3,134,  XN ( NO ) , EX (NO ) , YN ( NO ) , EY ( NO ) ,THK ( NO ) , 

1  ET(N0)  ,VN(NO) ,EV(NO) 

WRITE  OUTPUT  TAPE  3,135 

WRITE  OUTPUT  TAPE  3,136,  XO  ,Y0 , THDEG0, VO 

WRITE  OUTPUT  TAPE  3,137 

WRITE  OUTPUT  TAPE  3,139 

NOO  =  NO  -  I 

no  5059  KE  =2, NOO 

WRITE  OUTPUT  TAPE  3,140,  TN (KE ) ,XN ( KE ) , YN I KE ) ,THN ( KE ) , VN I KE ) , 

I  EX(KE),EY(KE),ET(KE),EV(KE) 

LINE  =  LINE  *  3 
IF  (56-LINE)  5058,5058,5059 

5058  WRITE  OUTPUT  TAPE  3,138 
WRITE  0UTPUT  TAPE  3,139 
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k 


LINe  *  5 

5059  C0NTINUE 

G0  T0  (5060t5060t5060t5060,5061)tN 

5060  WRITE  0UTPUT  TAPE  StlAl,  XLiN0 1  XO0S  i  XN  ( NI0)  t  YF 

5061  NJR  =  NJR  ♦  I 
JN  =  JN-l 

IF  (JN)  l.Lf2 
C 

END 
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APPENDIX  II 


This  appendix  contains  output  for  the  sample  calculation 
depicted  by  Figure  2.  The  input  for  this  computation  is  illustrated 
in  Section  C.  The  first  portion  of  the  output  consists  of  a  table  of 
range,  altitude,  speed  and  angle  as  a  function  of  time.  Also  in¬ 
cluded  are  the  time  derivatives  of  these  quantities,  the  mach  number, 
the  drag  coefficient,  the  thrust,  vehicle  weight  and  drag,  the  air 
density,  and  the  air  pressure  when  required. 

The  second  section  of  the  output  is  devoted  to  the  error  analysis 
routine.  This  begins  with  a  description  of  the  error  specification  and 
provides  the  estimated  error  for  comparison.  A  table  of  range  error 
(feet),  altitude  error  (feet),  speed  error  (feet/ second),  and  angular 
error  (degrees)  as  a  function  of  time  follow.  This  page  is  terminated 
by  bounds  for  the  impact  range  which  are  based  upon  a  consideration 
of  the  numerical  error  envelope  surroundin  g  the  computed  trajectory. 
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Copy  No. 

Feltman  Research  Laboratories 
Mr.  L.  H.  Eriksen,  Director  )  1 

Mr.  R.  Frye,  Deputy  Director  )  (in  turn) 

Engineering  Sciences  Laboratory 
Mr.  W.  R.  Benson,  Chief  2 

Computing  and  Analysis  Section 
Mr.  J.  Spector,  Chief  3-5 

Mathematical  Analysis  Unit 

Dr.  L.  F.  Nichols,  Chief  6-7 

Dr.  H.  J.  Kopp  8-17 

Mr.  B.  Barnett  18 

Digital  Applications  Unit 

Mr.  I.  Klugler,  Chief  19-28 

Mr.  J.  Nielsen  29  -  33 

Mrs.  E.  Stochel  34 

Digital  Systems  Unit 

Mr.  D.  L.  Grobstein,  Chief  35 

Special  Assistant  to  Chief,  CAS 
Dr.  S.  L.  Gerhard  36 

Aeroballistics  Section 

Mr.  A.  A.  Loeb,  Chief  37  -  38 

Analysis  Unit 

Mr.  D.  Mertz,  Chief  39  -  44 

Engineering  Research  Section 
Mr.  G.  Demitrack,  Chief  45 

Engineering  Analysis  Unit 

Dr.  F.K.  Soechting,  Chief  46 
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DISTRIBUTION  LIST  (Cont) 

Copy  No. 

Pyrotechnics  Laboratory 

Mr.  S.  Sage,  Chief  47  -  48 

Technical  Information  Section 

Mr.  M.  Costello,  Chief  49  -  54 

Nuclear  Weapons  Development  Division 
Mr.  M.  C.  Epton,  Chief  55 

Long  Range  Atomic  Warheads  Laboratory 
Mr.  H.  W.  Painter,  Chief  56 

Tactical  Atomic  Warheads  Laboratory 
Mr.  G.  I.  Jackman,  Chief  57 

Atomic  Ammunition  Development  Laboratory 
Mr.  J.  G.  Drake,  Chief  58 

Nuclear  Concepts  and  Systems  Section 
Mr.  M.  Rosenberg,  Chief  59 

Mr.  E.  Jaroszewski  60 

Mr.  F.  Scerbo  6l 

Ammunition  Development  Division 
Mr.  R.  Vogel,  Chief  62 

Warheads  and  Special  Projects  Laboratory 
Mr.  F.  E.  Saxe,  Chief  63 

Analytical  Section 

Mr.  S.  E.  Einbinder,  Chief  64 

Artillery  Ammunition  Laboratory 
Mr.  E.  H.  Buchanan,  Chief  65  -  66 

Solid  Rocket  Propulsion  Laboratory 
Mr.  A.  LoPresti,  Chief  67 

Mr.  S.  Wasserman,  Aeroballistics  Specialist  68  -  72 

Mr.  L.  Goldsmith  73 
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DISTRIBUTION  LIST  (Cont) 


Liquid  Rocket  Propulsion  Laboratory 
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